Hubbard-model description of the high-energy spin-weight distribution in La 2 Cu04 
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The spectral-weight distribution in recent neutron scattering experiments on the parent compound 
La2CuC>4 (LCO), which are limited in energy range to about 450 meV, is studied in the framework of 
the Hubbard model on the square lattice with effective nearest-neighbor transfer integral t and on- 
site repulsion U. Our study combines a number of numerical and theoretical approaches, including, 
in addition to standard treatments, density matrix renormalization group calculations for Hubbard 
cylinders and a suitable spinon approach for the spin excitations. The latter spin-1/2 spinons are the 
spins of the rotated electrons that singly occupy sites. These rotated electrons are mapped from the 
electrons by a uniquely defined unitary transformation, in which rotated-electron single and double 
occupancy are good quantum numbers for finite interaction values. Our results confirm that the 
U /8t magnitude suitable to LCO corresponds to intermediate U values smaller than the bandwidth 
8t, which we estimate to be 8t m 2.36 eV for U/8t w 0.76. This confirms the unsuitability of the 
conventional linear spin-wave theory. Our theoretical studies provide evidence for the occurrence 
of ground-state d-wave spinon pairing in the half-filled Hubbard model on the square lattice. This 
pairing applies only to the rotated-electron spin degrees of freedom, but it could play a role in a 
possible electron d-wave pairing formation upon hole doping. We find that the higher-energy spin 
spectral weight extends to about 566 meV and is located at and near the momentum [7r,7r]. The 
continuum weight energy- integrated intensity vanishes or is extremely small at momentum [it, 0]. 
This behavior of this intensity is consistent with that of the spin waves observed in recent high- 
energy neutron scattering experiments, which are damped at the momentum [n, 0]. We suggest that 
future LCO neutron scattering experiments scan the energies between 450 meV and 566 meV and 
momenta around [7r,7r]. 

PACS numbers: 78.70.Nx, 74.72.Cj, 71.10.Fd, 71.10.Hf 
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I. INTRODUCTION 

The development of a better understanding of quan- 
tum magnetism is important for improving our under- 
standing of the high-temperature cuprate superconduc- 
tors. Indeed, the parent compounds of the cuprates are 
insulating antiferromagnets, and these less complicated 
undoped systems can provide valuable information on 
which model Hamiltonians quantitatively describe the 
cuprates. Improved determination of the model Hamilto- 
nians is essential because of the many nearby competing 
phases in the doped systems, easily affected by small pa- 
rameters, which can now be seen because of continued 
improvements in numerical simulations 1 -. 

The spectral-weight distribution in recent neu- 
tron scattering experiments on the parent compound 
La2Cu04 (LCO), which are limited in energy range to 
about 450 meV, raise new interesting questions 2 -. In 
LCO, antiferromagnetic order occurs with a commensu- 
rate wave vector [n, n], where [it,tt] is observed to re- 
main commensurate for a finite level of doping. A [n, n] 
Goldstone mode was predicted by a spin-bag model 3 -. A 
decade ago the neutron scattering experiments on LCO 
of Coldea et alA first showed sufficient details of the spin- 
wave spectrum to demonstrate that a simple nearest- 
neighbor Heisenberg model must be supplemented by a 



number of additional terms, including ring exchanges. 
These terms arise naturally out of a single band Hub- 
bard model with finite U/t, and several detailed studies 
showed that the spin-wave data in the available energy 
window could be successfully described by the Hubbard 
model using a somewhat smaller value of U/t ~ 6 — 8 
than originally thought appropriated— . (For the effec- 
tive Coulomb repulsion U in units of the bandwidth, 8t, 
this refers to intermediate values, U/8t ~ 0.75 — 1.) 

Part of the spin spectral weight reported in Ref^ was 
deduced to be outside the energy window. The recent 
improved neutron scattering experiments of Refi 2 -, with 
a wider energy window of about 450 meV, have raised a 
number of questions. Surprisingly, these studies revealed 
that the high-energy spin waves are strongly damped 
near momentum [w, 0] and merge into a momentum- 
dependent continuum. These results led the authors of 
Refi 2 - to conclude that "the ground state of La2Cu04 con- 
tains additional correlations not captured by the Neel- 
SWT [spin-wave theory] picture". 

This raises the important question of whether the more 
detailed results can still be described in terms of a sim- 
ple Hubbard model. We show that the Hubbard model 
does describe the new neutron scattering results. Our 
results confirm that the U/t value suitable to LCO is in 
the range U/t £ (6,8). Inclusion of second- and third- 
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neighbor hopping parameters, t' and t" , into the Hubbard 
Hamiltonian lead to an interaction strength U/t w 8. 
Specifically, the studies of Ref4 have considered that the 
best fits to the ensemble of LCO inelastic neutron scat- 
tering points from Ref4 are reached for U/t rj 7.9 if one 
includes four independent parameters, t, t 1 , t", U, and for 
U/t « 7.1 if one includes only t and U. Furthermore, the 
results of Refs4^ reveal that as far as the LCO inelastic 
neutron scattering is concerned the inclusion of t 1 and t" 
does not lead to a better quantitative fit. Accordingly, 
the studies of this paper consider the half-filled Hubbard 
model on the square lattice with only two independent 
effective parameters, t and U. 

Our study uses a combination of a number of nu- 
merical and theoretical approaches, including, in ad- 
dition to standard treatments, density matrix renor- 
malization group (DMRG) calculations for Hubbard 
cylinders^— and, since conventional linear spin-wave 
theory is unsuitable, a spinon operator approach for the 
spin excitations^^. The spinon operator approach is 
suitable for LCO's intermediate U/t values, and corre- 
sponds to a particular case of a general operator rep- 
resentation that profits from the recently found model's 
extended global symmetry^. 

An exact result valid for the Hubbard model on any bi- 
partite lattice is that for onsite interaction U ^ it has 
two global SU (2) symmetric a 15 ' 16 , which refer to a global 
SO (4) = [SU(2) (g) SU(2)]/Z 2 symmetryil^. A recent 
study of the problem by one of us and collaborators re- 
ported in Ref J^, reveals that an exact extra global c hid- 
den U(l) symmetry emerges for U ^ 0, in addition to the 
50(4) symmetry. Specifically, the Hubbard model on a 
bipartite lattice, such as the present square lattice, has a 
global [SU{2)®SU(2)®U{l)]/Z 2 = [SO{4)®U(l)]/Z 2 = 
SO(3) <g) SO(3) ® 17(1) symmetry. The index c in the des- 
ignation c hidden symmetry is intended to distinguish it 
from the 77-spin U(l) symmetry and spin U(l) symme- 
try in the corresponding two model's SU(2) symmetries. 
The index c also labels the c fermions, whose occupancy 
configurations generate the representations of the global 
c hidden U{1) symmetry algebra. That the latter sym- 
metry is hidden follows from the fact that except in the 
U/t — > 00 limit, its generator does not commute with the 
electron - rotated-electron unitary operator. As a result, 
for finite U ft values its expression in terms of electron cre- 
ation and annihilation operators has an infinite number 
of terms. (The generators of the 77-spin and spin SU{2) 
symmetries commute with that unitary operator.) 

The origin of the extended global symmetry is a local 
gauge SU{2) <g) SU(2) <g> U(l) symmetry of the model 
Hamiltonian electron-interaction term first identified in 
RefJ^. That local symmetry becomes for finite U and 
t a group of permissible unitary transformations. The 
corresponding local U{1) canonical transformation is not 
the ordinary U(l) gauge subgroup of electromagnetism. 
It is rather a "nonlinear" transformation^. 

For very large U/t values the Hubbard model may be 
mapped onto a spin-only problem whose spins are those 



of the electrons that singly occupy sites. However, for in- 
termediate U/t values this mapping generates many com- 
plicated terms in the Hamiltonian, when written in terms 
of electron creation and annihilation operators. Here we 
address that problem by expressing the Hamiltonian in 
terms of the rotated-electron operators, which naturally 
emerge from the generators of the model's symmetries. 

In contrast to electrons, for rotated electrons single 
and double occupancy are good quantum numbers for 
U/t > 0. For large U/t values electrons and rotated 
electrons are the same objects. Apparently, the Hamil- 
tonian t/U expansion is formally similar in terms of elec- 
tron and rotated-electron operators. However that is 
only so for very large U/t values. For instance, there 
are well-defined t 2 /U and t 4 /U 3 terms in the Hamilto- 
nian expression in terms of rotated-electron operators, 
which are identical in form for very large U /t to the cor- 
responding terms using electron operators. However, for 
intermediate U/t <G (6,8), if one expresses the former 
t 2 /U and t 4 /U 3 terms in electron creation and annihi- 
lation operators, one finds many complicated higher or- 
der P /Ui~ l terms where for the half-filled case j are 
even integers, j = 2,4,6,... and j = 4,6,8,..., respec- 
tively. Hence the first few terms of the Hamiltonian ex- 
pression in terms of rotated-electron operators describe 
many higher-order electron processes. For moderate U/t 
the rotated-electron operators also generate a much sim- 
pler form for the energy eigenstates as well as for com- 
plicated processes involving a large number of electrons. 
Our spin- 1/2 spinons correspond to the spin- 1/2 spins 
of the rotated electrons that singly occupy sites, so that 
they are well defined for U/t > 0. 

When one decreases U/t to the intermediate U/t val- 
ues suitable for LCO, the above mentioned Hamilto- 
nian terms become increasingly important and generate 
higher-order spinon processes. Fortunately, those are 
simpler than the corresponding electron processes. In- 
deed, the use of our operational representation renders 
the intermediate U/t quantum problem in terms of ro- 
tated electrons similar to the corresponding large-t/ /t 
quantum problem in terms of electrons. The effect of 
decreasing U/t is mostly an increase of the energy band- 
width of an effective band associated with the spinon oc- 
cupancy configurations. The intermediate U/t rotated- 
electron processes may be associated with exchange con- 
stants describing rotated-electron motion touching pro- 
gressively larger number of sites. Within our Hubbard 
model's representation such Hamiltonian terms emerge 
naturally upon decreasing the magnitude of U ft. 

Our theoretical studies provide evidence of the occur- 
rence of ground-state d-wave spinon pairing in the half- 
filled Hubbard model on the square lattice. One of the 
few exact theorems that applies to the half-filled Hub- 
bard model on a bipartite lattice with a finite number 
of sites and thus on a square lattice is that its ground 
state is a spin-singlet stated. Within our spinon repre- 
sentation the ground-state spin-singlet A-spinon configu- 
ration corresponds to N/2 independent spin-neutral two- 
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spinon configurations. Under the spin-triplet excitation 
one of the N/2 spin-singlet spinon pairs is broken. Quan- 
titative agreement with the spin- wave spectrum obtained 
from our standard many-particle diagrammatic analysis 
is reached provided that the broken spinon pair has d- 
wave pairing in the initial ground state. Such a pairing 
refers only to the rotated-electrons spin degrees of free- 
dom. However, it could play a role in a possible d-wave 
electron pairing formation upon hole doping. 

Applying our approach to the new LCO high-energy 
neutron scattering reported in Refj£, we find that at mo- 
mentum [tt, 0] the continuum weight energy-integrated 
intensity vanishes or is extremely small. Furthermore, we 
find that beyond 450 meV, the spectral weight is mostly 
located around momentum [it, tt] and extends to about 
566 meV, suggesting directions for future experiments. 

The paper is organized as follows: In Sec. [II] the Hub- 
bard model on the square lattice and the basic quan- 
tities of our study are introduced. The description of 
the model's antiferromagnetic long-range order is ad- 
dressed in Sec. Mil Sec. IIVI presents a random-phase- 
approximation (RPA) study of the model's coherent spin- 
wave spectrum and intensity. Quantitative agreement 
with that observed in neutron-sattering experiments is 
used to find the U and t values suitable to LCO. The 
rotated-electron description emerging from the Hubbard 
model on the square lattice with extended global sym- 
metry is introduced in Sec. [V] This is the only section 
where general electronic densities and spin densities are 
considered. The goal of this more general analysis is the 
introduction of a spinon representation suitable to the 
LCO intermediate values of U/t. In Sec. IVII this spinon 
representation is used in the study of the general spin- 
triplet spectrum of the half-filled Hubbard model on the 
square lattice, which includes both the spin waves and 
the incoherent spin-weight continuum distribution. The 
comparison of the predicted spectral weights with those 
observed in the LCO high-energy neutron scattering is 
the goal of Sec. IVII1 Finally, Sec. IVIIII contains the con- 
cluding remarks. 



II. THE MODEL AND THE BASIC 
QUANTITIES OF OUR STUDY 

Most of our results refer to half-filling, so that the num- 
ber of lattice sites, N a , equals the number of electrons N. 
The exception is the general analysis reported in Sec. IVl 
which considers arbitrary values of the electronic density 
n = N/N a . The Hubbard model on a square lattice with 
N a 3> 1 sites and periodic boundary conditions reads, 

H = tf+UV D) 

f = = ~ Yl J2( C k,a + C \,,a <*i,°) > 



Vd = J2 (%>t - 1/2) (%4 - 1/2) • (1) 

3=1 



Here T is the kinetic-energy operator in units of t, Vd 
is the on-site repulsion interaction operator in units of 
U, ct and c?. a are electron creation and annihilation 
operators with site index j = 1, ...,N a and spin a =t,i, 
and npj.cr = ct CT Cf jiCT . The on-site repulsion interac- 
tion operator Vd may alternatively be expressed in terms 
of the electron double-occupancy operator D or single- 
occupancy operator Q given by, 

D = {N-Q)/2 ; Q = J2I2 ^-H,-*) > ( 2 ) 

3=1 a=f,l 

respectively. The expectation values, 



3 = 1 
N, 



(l-2d) = ^-f^iGSKh^-h^flGS), 



m AF 



3 = 1 

1 Na 1 

JVa 3 = 1 



[1 - 2SS] m° AF , 



(3) 



play an important role in our study, following the strong 
evidence that for U > and N a —> oo the model's ground 
state has antiferromagnetic long-range order—. In the 
last expression of Eq. ([3]) one has that j is an even inte- 
ger and an odd integer for each of the two sub-lattices, 
respectively. Moreover, in that expression m° AF stands 
for a mean-field sub-lattice magnetization that does not 
account for the effect of transverse fluctuations while SS 
does account for this effect, its value being estimated be- 
low. Specifically, m° AF is the sub-lattice magnetization 
of the spin-density wave (SDW) state obtained in a stan- 
dard mean-field treatment of the Hubbard interaction at 
zero absolute temperature as given, for instance, in Fig. 
3 of Ref£. 

The on-site spin operators involved in our studies read, 



r 3 ,s 2i l r i' 



1, 



(4) 



S r 3 ,s = —j[ n ^^~ n ^i 

The index s in these operators distinguishes them from 
the corresponding local operators associated with the rj- 
spin SU (2) symmetry algebra considered below in Sec. 

m 

Our study also involves the spin dynamical structure 
factors, 

V ; N a .4- 

3:3 =1 

j* OO 

x / dte™{GS\% jia (t)ig (0)\GS), (5) 



4 



Analytical approximation 

t/U expansion 

i DMRG, 10x6 cylinder 
o DMRG, 10x4 cylinder 



U/t 



Figure 1: Average single-occupancy: approximate expression 
[1 +tanh([//8t)]/2 valid for U/t < 8 (solid line), the limiting 
[//( > 1 expression [1 — cc(8t/U) 2 ] (dashed line), and from 
DMRG numerical results on two different width cylinders. 



where 



= x, y, z or a = — and a' = + and be- 



low we consider g = 2. It is straightforward to show 
that the sum rules [1/N a ] J2k S aa ' (k) where S aa ' (k) = 
[1 /2tt] dw S aa (jfe, tu) involve the average single occu- 
pancy (1 — 2d) and read, 



-E 

k 



S aa (k) 



[5 ata > + 2S a ,-Sa',+]{l ~ 2d) 



(6) 

In an ideal experiment all components are detected 
with equal sensitivity. As discussed for instance in Ref£, 
in that case a transfer of spectral weight from the lon- 
gitudinal to the transverse part as the energy increases 
is observed. Hence independent of the scattering geome- 
try, the corresponding effective spin dynamical structure 
factor satisfies the sum rule, 



1 



1 



k 



— > - du S e ^(k, u) = /4 2(1 - 2d) . (7) 



That the coefficient involved is 2(1 — 2d) rather than 3(1 — 
2d) follows from one mode being perpendicular to the 
plane and thus silent in the experiment^. 



III. THE HUBBARD MODEL ON THE SQUARE 
LATTICE ANTIFERROMAGNETIC 
LONG-RANGE ORDER 

For the range U/t G (0, 8), the antiferromagnetic long- 
range order may be accounted for by a variational ground 
state with a SDW initial trial state, such as for instance 
a Gutzwiller projected antiferromagnetic stated, 

\G) = e- gb \SDW) , U/t<8, (8) 
or the following related state, 

\GB) = e- ht ^e- st> \SDW) , U/t < 8. (9) 



Here \SDW) is the ground state of a simple effective 
mean-field Hamiltonian, such as that of Eq. (18) of 
Refill. For > 1 this order is as well accounted for 
by a Baeriswyl variational state, 



|B> = e- h# / t |c»>, Etyt»l 



(10) 



where |oo) is the exact U/t — s- oo ground stated. The 
coefficients h and g multiplying the kinetic-energy and 
double-occupancy operators, respectively, in the state ex- 
pressions given in Eqs. (j8j)- (fTU)) are variational parame- 
ters. The above states involve as well a variational gap 
parameter A, which is expected to tend to zero as the 
trial state approaches the exact ground state. Indeed, 
that variational parameter is an infinitesimal symmetry- 
breaking field. 

Similarly for the trial state \SDW), the relation 



a Waf? = (! - 4d ) 



(11) 



holds for the states |G), \GB), and \B). However, the 
corresponding function d = d(U/t) is in general state 
dependent. Inversion of the simple relation provided in 
Eq. (|XTJ) gives, 



-[1 - 4[m 



AFl 



(12) 



This is consistent with d not being affected by transverse 
fluctuations. 

The evaluation of the ground-state energy for \G) and 
\GB) is for N a > 1 an involved problem. Here we resort 
to an approximation, which corresponds to the simplest 
expression of the general form, 



E/N = T oqu + Ud; 



(13) 



compatible with three requirements. Those are: 



1) The relation d = 
12) must be fulfilled: 



|[l-4[m 



AFl 



provided in Eq. 



2) The antiferromagnetic long-range order must occur 
for the whole U/t > range; 

3) The small- U ft expansion of the energy E/N of Eq. 
(|13p must lack of a linear kinetic-energy term in U for 
U/t <C 1 (except for the term U d corresponding to the 
on-site repulsion.) 

Brinkman and Rice found qu = 8d(l — 2d) for the 
original paramagnetic-state Gutzwiller approximation^, 
which is lattice insensitive and thus does not account 
for the square-lattice antiferromagnetic long-range order. 
The simplest modified form of the quantity qu suitable 
to a broken-symmetry ground state such that the above 
three conditions are met is, 



a\ +) d 



(1 - 2d) 
4[m°, 



a-2 



AFl 



a 3 ■ 



(14) 
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Here. 



a[ ±} =tt 2 ±4; 



o a = (1 - br 2 /2a[ +) ]) ; 



(15) 



and the coefficient 03 is a function 0,3 = a 3 (U/t) of Z7/t 
whose approximate limiting behaviors are, 



"3 



,(-) r 



U, 



.anh^p + ai+^M" 3 
8t 



= -coK/2] 2 -, C//t>l 



where, 



Co 



a 1 

4 + 8 



0.1462. 



U/t<8, 
(16) 

(17) 



states |G) and \GB) (see Fig. 4 of Ref£). For > 1 
we find the behavior 4[m^] 2 w e- 2c <> for the state 

\B), so that, 



(1 - 2d) 



1 + tanh 



, U/t < 



1 - c 



, U/t > 1 . 



(19) 



Furthermore, the states ISDVF) and |G) give a sub- 
lattice magnetization ttiaf ~ "Af = 5VI — 4d, with 
an improved U/t dependence, d w j[l — tanh(£//8t)], 
for the latter, as follows from the corresponding (1 — 2d) 
expression of Eq. (|T9"|) . For the state \G) the sub- lattice 
magnetization is then given by, 



The corresponding estimate a = 0.6696 is that of the 
Heisenberg-model studies of Refj 2 ^. Moreover, the quan- 
tity 4[m AF } 2 on the right-hand side of Eq. (|T4"|) behaves 



L AF 




tanh 



U/t < 



(20) 





AFl 



U/8t for U/t<l. 



as 4[m 

Minimization of the ground-state energy defined by 
Eqs. (fT5)) - (fl"B]) with respect to d leads indeed to d = 



On the other hand, we find that the states \GB) and 
\B) have m AF and rn AF sub-lattice magnetization nu- 
merical values very close to those given by the relation 



\[l - 4[m° F ] 2 ]- The limiting behaviors of that energy ^ _ 2SS] m ^ of Eq p ^ 



4 

are 



E/N « To + Ud 
8f 2 



1 [/ 2 

8^T' ^ <<X 



4c ^, Ltyt»l 



(18) 



We note that the small-U/t second-order coefficient reads 
-[1/8tt 2 ] « -0.0127, in agreement with that, « -0.0127, 
obtained by second-order perturbation theory^ 3 .. For 
U/t ^ 1 one recovers the known result E/N = 
— 4co[8t 2 /?7} 2 ^, so that our approximation agrees with 
the known limiting behaviors. 

The quantity 4[m^ F ] 2 on the right-hand side of Eq. 
([Till behaves as 4[m AF ] 2 = U/8t for U/t < 1. However, 
its U/t dependence for the range U/t € (0, 8) remains an 
open problem. Here we have performed DMRG calcula- 
tions of (1 — 2d), which according to the relation provided 
inEq. (12) is given by (1-2(2) = ±(1 + A[m AF ] 2 ). Hence 
its U/t dependence fully determines that of 4[m^ F ] 2 . 
The corresponding DMRG results are shown in Fig. 1. 
Specifically, two different circumference cylinders were 
simulated as a function of U/t, with open boundary con- 
ditions in x and periodic in y, and the double occupancy 
measured in one of the middle columns. A maximum of 
in = 4000 states were kept, with an accuracy of ~ 10~ 4 
in (1 — 2d) for the 10 x 4 system for the least accurate 
smaller U /t values, and about 10 -3 for the 10 x 6 system. 
We find that the value of (1 — 2d) is relatively insensitive 
to cluster size, and these cluster sizes are representative 
of two-dimensional (2D) behavior—. 

We then find that 4[m AF ] 2 ps tanh(C//8t) gives for the 
range U/t G (0,8) quantitative agreement for the (1 — 
2d) dependence on U /t with both our numerical DMRG 
calculations (see Fig. 1) and the numerical results for the 



SS 



d, U/t<8, 



m H AF 



m 





HAF 



c//t> 1 



(21) 



respectively. Here TTfj 1AF = 1/2 and tuhaf r 
0.303 is the Heisenberg-model's sub-lattice magnetiza 
tion magnitude 2 ^, so that 6S w d + 0.197 in Eq. (J2T) 
for U/t^> 1. Hence one finds, 



L AF 



rn B AF 



1 + tanh — 
8t 



0.303 - 0.803 x c ( — 




( ^ ) 



E//*>1. (22) 



The magnitudes of the sub-lattice magnetizations m AF 
and m AF as given in Eqs. (|2U)l and (|2"2"f for the states 
\G) and \GB), respectively, are provided in Table Q] for 
several U /t values. In that table the magnitudes of a 
sub-lattice magnetization m l AF er that for U/t > we de- 
fine as m l AF er = (1 — 2d) [rriHAF /^haf} 171 af are a ^ so 
given. Note that for U/t 3> 1 the sub-lattice magnetiza- 
tion m l AF er becomes m AF . Probably it is closest to the 
exact mAF, while m AF is that consistent with our use of 
the RPA in the ensuing section, to study the spin-wave 
spectrum and corresponding intensity. 



IV. COHERENT SPIN- WAVE SPECTRUM AND 
INTENSITY AND LCO U AND t VALUES 

To study the coherent spin-wave weight distribution 
and spectrum, we have calculated the transverse dynam- 
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u/t 


6.1 


6.5 


8.0 


10.0 


m AF 


0.401 


0.410 


0.436 (0.43&2i) 


0.461 (0.456^) 


GB 

m AF 


0.329 


0.342 


0.384 (0.39^.) 


0.426 


lower 

m AF 


0.200 


0.207 


0.233 


0.258 



Table I: The sub-lattice magnetizations as calculated here for 
several U/t values and some results from Refs.— 



ical susceptibility, 



X - + (k,r) 



E e -*(fj-^) 

... 

x (s~Jr)s+ r (0)), (23) 



in the RPA. Here r denotes the imaginary time in Mat- 
subara formalism and we shall take the zero temperature 
limit. Because we deal with the antiferromagnetic order 
(Neel state) it is convenient to define two sub-lattices, a 
and b, and redefine the susceptibility as a 2x2 tensor x^,v 
where the Greek subscripts denote sub-lattice indices, 



(Na/2) 



E 



W.,,.(o)> 



(24) 



The original susceptibility in Eq. (|23[) is then simply 
related to this tensor as, 

X(k,T) = ^[Xaa+Xbb + Xab+Xba] ■ (25) 

We define electron field operators for each sub-lattice, 
fir and br , as, 

KG KG 1 1 



Crj(Eb,(T 



-7^= r 



1 



E 



ko 



ko 



(26) 
(27) 



In the momentum summations the reduced Brillouin zone 
(RBZ) covers only half of the original Brillouin zone (BZ) 
for the square lattice. The effective Hamiltonian that 
describes the SDW phase in mean field theory for the 
Hubbard interaction can be written as, 



H eff 

with, 



yffit hi 

\ ka ko 



^ f(k) 
f(k) -e* 




(28) 



e<r = -U — m AF , 



f(k) 



-2t [cos(k x ) + cos(fc y )] 



(29) 



Using the effective Hamiltonian given in Eq. (|28[) one 

derives from Eq. (|24p a susceptibility tensor x^i(fc, r) 
at mean field theory level. 



Treating the Hubbard interaction futher in the RPA 
and Fourier transforming the susceptibility from imag- 
inary time t to (k,ioj) space, the susceptibility tensor 
then obeys the Dyson equation, 



-RPA 
A 



x (0) 



Ux {0) x RP 



which can be recast as, 

= [f-UxW]- 1 *™ 



-RPA 



(30) 



(31) 



Here / stands for the 2x2 identity matrix. Such a pro- 
cedure of treating the interaction in RPA on top of the 
mean field solution has been used in previous studies^^. 

X RPA has a pole ioj = u>(k) obtained from the equation 
Det [1— Ux ] = 0, which provides the dispersion relation 
uj(k) for the spin waves. It has been shown in Ref£ that 
an excellent agreement with the spin- wave spectrum from 
Refi^ is achieved. In Fig. [5] (top) we show a fit to the 
more recent experimental data of Ref (solid line) along 
with the results from the si fermion method reported 
below in Sec. IVII (dashed line) for U/t = 6.1 and t = 
295 meV. This corresponds to a bandwidth 8t rj 2.36 eV. 

Importantly, provided that the t magnitude is slightly 
increased for increasing values of U /t, agreement with the 
LCO spin-weight spectrum and distribution can be ob- 
tained for the range U/t G (6, 8) and thus U/8t G (0.75, 1) 
in units of the bandwidth 8t. For U /t values smaller than 
6 (and larger than 8), the spin- wave dispersion between 
[it, 0] and [tt/2, 7r/2] has a too large energy bandwidth 
(and is too fiat) for any reasonable value of t. 

Let \u, cj(fc)) denote the excited energy eigenstates of 
energy uj(k) and momentum k that contribute to the 
coherent spin-wave spectral weight. In the case of the 
— h spin dynamical structure factor given in Eq. ([5]) for 
a = — and a' = +, the corresponding coherent spin- wave 
spectral weight in units of ijl 2 b is given by Zd2(l — 2d). 
The factor Zd in this expression reads, 



1 



N a (l 



2d) E E i 



(32) 



where st is the Fourier transform of the spin operator 

k,s 

si s defined in Eq. ^ and the sum over energy eigen- 
states excludes those that generate the coherent spin- 
wave weight, \v) = \v, uj{k)). In the U/t — > oo limit, Zd 
may be identified with the corresponding Zd = Z c Z x 
factor of the Heisenberg model on the square lattice. 
According to the results of Ref.—, the factors Z c and 
Z x have magnitudes Z c w 1.18 and Z x w 0.48, respec- 
tively, so that Zd ~ 0.57. The limiting values Zd = 1 for 
U/t — > and Zd ~ 0.57 for U/t — > oc and the approx- 
imate intermediate valued, Zd ~ 0.65, at U/t = 8 are 
recovered as solutions of the equation, 



„ -Z d tanhl V 4TT 

Z d = e V v 4xt t 

which is used here for finite U /t. 



(33) 
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Figure 2: (Top) Spin-wave excitation spectrum along BZ spe- 
cial directions as specified in Ref.—. (Bottom) Spin-wave in- 
tensity as obtained from the poles of the susceptybility (see 
text). Experimental points from Ref.—. 



The spin dynamical structure factor measured in the 
high-energy inelastic neutron scattering experiments of 
Refj 2 - includes a Bragg peak associated with its elastic 
part. The corresponding elastic spectral weight is in- 
cluded in the total spin- weight sum-rule, y? B 2(1 — 2d), 
of Eq. ([7]). In the thermodynamic limit the upper- 
Hubbard band processes generate nearly no spin weight. 
Hence the longitudinal spectral weight within the sum- 
rule fi B 2(1 — 2d) refers to the elastic contribution. The 
elastic weight is given by ps y? B 4(to j 4f) 2 - The inelastic 
spin spectral weight corresponds to the remaining weight 
in the spin-weight sum-rule ps fi B [2(1 — 2d) — A(mAF) 2 ]- 
It refers to the — h spin dynamical structure factor given 
in Eq. for a = — and a' = +. It then follows that 
the experimentally determined spin- wave intensity, which 
corresponds to the coherent part of the inelastic spin 
spectral weight, is in units of fj, B approximately given 
by, 

W SW « Z d [2(1 - 2d) - 4(m AF ) 2 ] . (34) 

The GA+RPA method used in Ref£ accounts for the 
quantum fluctuations that control the longitudinal and 
transverse relative weights. Within our description no- 
tations, that method is designed to make the inelastic 
spin spectral weight fi B [2(1 — 2d) — 4(m J 4 F ) 2 ] rather than 
fx% 2(1 - 2d). Hence for that GA+RPA method the spin- 
wave intensity factor is Zd as defined in Eq. (|32l) . 

On the other hand, the RPA used here refers to the — h 
spin dynamical structure factor alone. Hence it implic- 



itly considers that the total inelastic spin spectral weight 
is Mb 2(1 - 2d) rather than /j| [2(1 - 2d) - 4(m j4F ) 2 ]. 
Therefore, to describe the actual spin- wave intensity mo- 
mentum distribution one must use a corresponding ex- 
perimentally determined factor Z d xp < Zd such that, 

W SW = Z e d xp 2(l-2d) = Z d [2{l-2d)-A(m AF f] . (35) 

From the residue of the spin- wave pole the susceptibil- 
ity coherent part then reads, 

-+tu ■ \ 7 ex P Res [x~+(fc, Zw(fc)] 

Xcc , = Z d " ^ ■ — jr , (36) 

;=±1 tuj — I ujyk) 

with x ^ obtained in RPA above. The measured inten- 
sity is21, 

I sw (k) = n[S xx (k) + S yv {k)] = nS~ + {k) . (37) 

In Fig. [5] (bottom), we plot the corresponding RPA 
spin-wave intensity, 

Isw(k) = -br/2]Zf p Res [*-+(£, a>(*)] . (38) 

The good agreement with the experimental data, spe- 
cially near the point M, reproduces the theoretical results 
of Ref? 2 .. It is here obtained for the value Z e d xp sa 0.49, 
which corresponds to the choice tuaf — m AF = m AF^~ 
2d) such that SS ps d. The m^p dependence on U jt is 
given in Eq. (J22J) for U/t < 8. As in Ref<2, the spin- 
wave intensity shows disagreement around the X point, 
which here probably stems from effects not captured by 
the RPA. 



V. THE ROTATED-ELECTRON DESCRIPTION 
EMERGING FROM THE MODEL'S EXTENDED 
GLOBAL SYMMETRY 

The goal of this section is the introduction of the 
general rotated-electron representation from which the 
spinon representation used in the ensuing section natu- 
rally emerges. In contrast to the remaining sections of 
this paper, here we consider arbitrary values of the elec- 
tronic density n = N/N a and spin density m = [Nf — iVj,] . 

A. The electron - rotated-electron unitary operator 

We denote the spin and 77-spin of an energy eigen- 
state by S s and S v , respectively. The corresponding 
spin and 77-spin projections read S* = — \ [iV-j- — AjJ and 
S* = — 5 [N a — N] , respectively. The lowest-weight states 
(LWSs) of both the 77-spin and spin algebras are such that 
S a = —S^ where a = 77 for 77-spin and a = s for spin. 
The numbers, 

n v = S v -^(N a -N) = 0,l,...,2S v , 

n s = S s -^(N t -N i ) = 0,l,...,2S s , (39) 
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vanish for such a LWS. 

Let {\^i rt i, ,.u)} be a complete set of A Na energy, 
momentum, ?7-spin, 7/-spin projection, spin, and spin- 
projection eigenstates for u = U/t > 0. Here l r]s is 
a short notation for the set of four quantum numbers 
[5*,,, 5* s , n r; , n s ] and the index l r represents all remaining 
quantum numbers, other than those, that are needed to 
fully specify an energy eigenstate \^i r .i v , lU ) ■ The energy 
eigenstates of that set that are not LWSs are generated 
from those as follows, 



I* 



l r ,lris ,U/ 



n 



i 



' Co. 



{sir- 



«) ■ (40) 



Here. 



C a = <*M° s ,«l(^) nQ (^)""l*M« s ,,) 

= [n a !] JJ [ 25 Q + 1 — j'] , a = V ,s, (41) 

for n a = l,...,2S a are normalization constants, the 77- 
spin (a = 77) and spin (a = s) off-diagonal generators 
and S a are given in Eq. (|A7[) of Appendix El and l ns and 
Z° s stand for [S r] , S s ,n v , n s ] and [S r] , S s ,0, 0], respectively. 
Within our notation, l° s refers to values of the general 
index associated with a LWS such that n J( = n s = 0. 

For the Hubbard model on the square lattice and also 
on the ID lattice, upon adiabatically increasing U ft from 
any finite value to the U/t — > 00 limit, each energy 
eigenstate |\&j ri j u ) continuously evolves into a uniquely 
defined corresponding energy eigenstate l 1 ^;,. 00)1 and 
vice versa. We emphasize though that due to the high 
degeneracy among different spin sectors as well as 77-spin 
sectors that occurs in the U/t — > 00 limit, there are in 
such a limit many more choices of energy eigenstates sets 
than for U jt finite. Accordingly, upon adiabatically de- 
creasing U jt most of such U/t 00 states do not evolve 
into finite-[/ /t energy eigenstates. Our above procedure 
uniquely defines a convenient set of U/t — > 00 energy 
eigenstates that upon adiabatically decreasing U/t do 
evolve into finite-?/ /t energy eigenstates. 

Both the corresponding sets of 4^" states {|'I , ;,.,z ? , s ,ti)} 
and {| v I / i r .z, !S ,oo)}, respectively, are complete and refer to 
the same Hilbert space. Hence there is a uniquely defined 
unitary transformation connecting the states \^i r ,i va ,u) 
and | x f , i ri i 7Js . 00 )- Indeed, since the model's Hilbert space is 
the same for all U/t > values considered here, it follows 
from basic quantum-mechanics Hilbert-space and opera- 
tor properties that for this choice there exists exactly one 
unitary operator V = V(U/t) such that any U/t — > 00 
energy eigenstate \'^i rt i ri3 .oc) is transformed onto the cor- 
responding U/t > energy eigenstate \^i rt i, B . u ) as, 



I* 



(42) 



The energy eigenstates |*i r ,i, s , u ) = \^i r ,i vs ,oo) (one 
for each value of U/t > 0) that are generated from the 



same initial U/t —> 00 energy eingenstate \^i ri i ai00 ) be- 
long to the same V tower. 

The rotated-electron operators are given by, 



ct 



(43) 



For U/t > the operator S appearing here can be ex- 
panded in a series of t/U whose leading-order term is 
provided in Eq. (|A2[) of Appendix |A"1 

Since the electron - rotated-electron unitary operator 
V commutes with itself, the equalities V = e~ s = V = 
e~ s and S = S hold. Hence both the operators V and 
S have the same expression in terms of electron and 
rotated-electron creation and annihilation operators. It 
then follows from the expression of the operator S pro- 
vided in Eq. (|A2[) of Appendix [A] that the corresponding 
rotated operator S has the following leading-order term, 



s = 'u 



2+1 - T_i 



(44) 



The rotated kinetic operators T +1 and T_\ appearing 
here and the related rotated kinetic operator T are given 
in Eq. (|A3[) of Appendix [S] The expressions of the cor- 
responding unrotated kinetic operators To, T+i, and T_i 
are provided in Eq. (|Al[) of that Appendix. 

Note that the equality S = S refers to the whole ex- 
pression of these operators. An important property for 
our study is that except in the U/t — > 00 limit the lead- 
ing order terms [f +1 - T_i] and [f +1 - T_i] of 
the operators S and S given in Eq. (|A2[) of Appendix 
lAl and Eq. (|44| . respectively, are different operators. 
Moreover, except for U/t — > 00 one has that To 7^ To, 
T + i 7^ and T_i 7^ T_i. This is behind for interme- 
diate U/t values the few first terms of the Hamiltonian 
t/U expansion as written in terms of rotated-electron op- 
erators containing much more complicated higher-order 
terms when expressed in terms of electron creation and 
annihilation operators. 

The main point here is that for the rotated electrons 
that emerge from the unitary transformation of Eq. (|4"5)) 
single and double occupancy are good quantum numbers 
for U/t > 0. Indeed, on any bipartite lattice the number 
of rotated-electron singly occupied sites operator, 



2S C = ViQv = J2s?i 



Sr^c = V 4 Sf^cV = ^2 fijt^cr (1 - -„) , (45) 
<r=T4 

commutes with the Hubbard model Hamiltonian^. Here 
Q is the corresponding number of electron singly occu- 
pied sites operator given in Eq. ©. This follows in 
part from the symmetries of the Hamiltonian electron- 
interaction term, which imply that all U/t — > 00 energy 
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eigenstates of the set {1^;,.,; s ,oo)} are as well eigenstates 
of the electron double-occupancy operator D and single- 
occupancy operator Q provided in that equation. Hence 
in the U/t — >• oo limit the Hilbert space is classified 
in subspaces with different numbers of doubly-occupied 
sites and each of the states {|^i r ,{_ s ,oo)} is contained in 
only one of these subspaces. The same applies to the 4 Na 
energy eigenstates of the set {|^'z r .,i, )5 ,«)} for U/t > in 
terms of rotated-electron doubly-occupied sites. 

The unitary operator V of our formulation is 
uniquely defined by its 4^" x 4 Na matrix elements, 
(*i r ,V,«|^|*i;,i^, u ). For U/t > most of these ma- 
trix elements vanish. For U/t — > oo rotated electrons 
become electrons so that the matrix representing the 
unitary operator V becomes the A Na x 4^" unit ma- 
trix. Hence (tf/^.oo |^|*z;^ s ,oo) = tfjpJ' 0n 
the other hand, as justified in Appendix [AJ the unitary 
operator V = V(U/t) commutes with the six generators 
of the global 77-spin and spin SU(2) symmetries. This 
implies that the matrix elements between energy eigen- 
states with different values of S^, S s , n v , and n s and thus 
of l vs vanish. Hence the finite matrix elements are be- 
tween states with the same l vs values so that we denote 
them by Vj r ,j/ , 

(*l r ,l n ,,oc\V\*V r ,l' n ,,°o) = K^V 1t , k , (46) 

where, 

= (^i r .i vs ,u\^v r .i vs ,oo) = (^ii.,i v ., 00 \%r,h.,u)* ■ (47) 

Given a complete set of 4^= energy, momentum, 77- 
spin, 77-spin projection, spin, and spin-projection eigen- 
states, {|^'; ri ; s .u)}, the electron - rotated-electron uni- 
tary operator considered here is for U/t > uniquely de- 
fined by the matrix elements of Eqs. (|46| and (|47|) . This 
corresponds to one out of the infinite choices of electron - 
rotated-electron unitary operators^. All these operators 
and corresponding unitary transformations refer to the 
same subspaces with fixed numbers of doubly-occupied 
sites, 0,1,2,3,.... They only differ in the choice of basis 
states within each of these subspaces. For most of these 
unitary operators the states V^\^i ri i VSj0 o) are not energy 
eigenstates for finite U/t values. The electron - rotated- 
electron unitary transformation considered here has been 
constructed to make these states energy eigenstates for 
finite U/t values, as given in Eq. (|42|) . 

B. The general operational description naturally 
emerging from the rotated electrons and symmetry 

The electron - rotated-electron unitary transforma- 
tion is closely related to the extended global S*0(3) x 
50(3) x U(l) symmetry found in RefJ^ for the Hamil- 
tonian given in Eq. ([1} on any bipartite lattice. Un- 
til recently^ it was believed that the model's global 



symmetry was for finite on-site interaction values only 
SO{4) = [SU{2)®SU{2)]/Z 2 . The occurrence of a global 
c hidden £7(1) symmetry beyond SO(A) in the model's 
global SO{3) <8> SO (3) <8 U(l) = [SO (4) ® U{l)]/Z 2 sym- 
metry must be accounted for in studies of the Hubbard 
model on any bipartite lattice. Such a global symme- 
try may be rewritten as [St/ (2) x SU{2) x t/(l)]/Z| and 
stems from the U ^0 local gauge £17(2) x SU(2) x t/(l) 
symmetry of the Hubbard model on a bipartite lattice 
with vanishing transfer integral, t = 0^. The seven local 
generators of the corresponding two gauge SU(2) sym- 
metries and t/(l) symmetry are the three spin local op- 
erators Sp. s provided in Eq. (j4]) and the three 77-spin 

local operators sj?. and the local operator sp^ c given in 
Eqs. (|A8j) and ([All|) of Appendix [A] respectively. The 
index I in the generators of the two SU(2) symmetries 
stand for I = ±, z. 

An important point is that although addition of 
chemical-potential and magnetic-field operator terms to 
the Hubbard model on a square lattice Hamiltonian 
given in Eq. ([T]) lowers its symmetry, these terms com- 
mute with it. Therefore, the global symmetry of the 
latter Hamiltonian being SO (3) <E) SO(3) <g> t/(l) im- 
plies that the set of independent rotated-electron occu- 
pancy configurations that generate all 4 Na energy eigen- 
states, {\^i r ,i ns , u }}, generate as well representations of 
the global symmetry algebra for all values of electronic 
density n and spin density m. It is confirmed in RefJ^ 
that the number of these independent [SU (2) <8> SU(2) ® 
U(l)]/Z$ = SO(3)®SO(3)®U(l) symmetry algebra rep- 
resentations equals for the present model on a bipartite 
lattice its Hilbert-space dimension, 4 Na . 

The generator sp. c of the local gauge U{1) symmetry 
given in Eq. (|A9j) of Appendix [A] and the alternative 
local generator c = (1 — s? jjC ) may be expressed as, 

Srj,c = Sfj = ff J>c frj,c j 

4„c = 0--& i )=M,cfi„e- ( 4 8) 

Here ft and f?. c stand for the following creation and 
annihilation operators, respectively, of suitable spin- less 
and 77-spin-less fermions, 

where we used that e" r ' rj = e~ l7r ' Tj . Here and throughout 
this paper the vector 7? has Cartesian components 7? = 

[7T,7r]. 

We call c fermions the rotated-electron related ob- 
jects whose creation and annihilation operators ft = 

^ ff } ,c V and fa," = ^ j^.c ^' respectively, are gener- 
ated from those of the spin-less and 77-spin-less fermions 
of Eq. (|4"9")l by the specific electron - rotated-electron 
unitary transformation uniquely defined by the matrix 
elements of Eqs. (|4"6")l and (|4"T1) . (No upper index / is 
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used onto the (rotated) c fermion operator f? )C .) Hence 
these operators read, 

4 lC = 4 3 , t t 1 - J + e "" j 5 ^ ,t ^ 4 . 

/^■,c = c^.tC 1 -«fv,4.) +e iff ' ?J ' ct t n ?3 4. (50) 

The rotated-electron creation and annihilation operators 
appearing here are generated from corresponding elec- 
tron operators by the unitary transformation uniquely 
defined above, as given in Eq. (|43| . The corresponding 
c fermion local density operator is given by, 

Qh=ficM,c (51) 

The c fermions live on a lattice identical to the original 
lattice. One can introduce c fermion momentum depen- 
dent operatorsZ^, 

4c=7^E e+i ^'4.c; ; =l >->^- ( 52 ) 

Here the c fermion operators /I where the index 

j' = 1, iV a refers to the sites of the original lattice are 
mapped from the rotated-electron operators by an exact 
local transformation given in Eq. f|50[) . The c momen- 
tum band has N a discrete momentum values qj where 
j = 1, N a . It has the same shape and momentum area 
as the electronic first-BZ. 

The generator 2S C of the related global c hidden U(l) 
symmetry in [SU(2) x SU(2) x U(l)]/Z$ found in Ref^ 
is the number of rotated-electron singly occupied sites 
given in Eq. (|45|) . Hence it involves the site summation 
y over the rotated local generator sp jjC rather than 
over the corresponding unrotated local operator s? jlC of 
Eq. (gU) of Appendix El This is why 2S C = V^QV, as 
given in Eq. (|45p . where the operator Q is that of Eq. 
@. The eigenvalues 2S C = 0, 1, 2, ... of the generator 2S C 
are thus the numbers of rotated-electron singly occupied 
sites. 

The c fermion creation and annihilation operators are 
found in Appendix [X] to obey the anti-commutation rela- 
tions given in Eq. (|A10|) of that Appendix. A straightfor- 
ward operator algebra then confirms that the c fermion 
local density operator of Eq. (fSTj) is the local operator 
Sf jjC appearing in the expression provided in Eq. (|45[) . 
Hence the global c hidden U(l) symmetry generator may 
be simply rewritten as, 

N a N a 
3=1 3=1 

One finds that except in the U/t —> oo limit the 
inequality J2j=i 9f ^ <2f holds. This confirms 

that the generator 2S C given in Eq. (|5"3")l does not com- 
mute with the electron - rotated-electron unitary opera- 
tor V = V. On the other hand and as justified in Ap- 



pendix [A] the three components of the momentum oper- 
ator P, three generators of the global spin SU{2) sym- 
metry, and three generators of the global 77-spin SU{2) 
symmetry commute with that unitary operator. Hence 
in contrast to the Hamiltonian and generator 2S Cl these 
operators have the same expression in terms of electron 
and rotated-electron creation and annihilation operators, 
as given in Eqs. ([A"6j) and (TAT)) of Appendix [A] On the 
contrary, the generator of the global c hidden U(l) sym- 
metry given in Eq. (|45[) does not commute with the uni- 
tary operator V. This is behind the hidden character of 
such a symmetry. 

Site summation Y]j—i over the rotated local operator 
provided in Eq. (j5"Tj) and over the following six rotated 
local operators, 







si 










i r _ - , 


°rj,s 




si 


= ,1 ,t ) 


sZ 

Tj.S 





gives the seven generators of the model's global 50(3) <E) 
50(3) ® 17(1) = [SU(2) ® SU(2) ® U(l)]/Z$ symmetry, 
as provided in Eqs. (j45|) and Eq. (|A7|) of Appendix |A"1 
However, except in the U/t — > 00 limit the six rotated 
local operators given in Eq. (|54p and the corresponding 
six unrotated local operators provided in Eq. ^ and Eq. 
(jA8[) of Appendix lAl are different operators. 

Interestingly, the 77-spin and spin SU (2) symmetries 
are within the present representation particular cases of 
a general r/s quasi-spin SU(2) symmetry. The corre- 
sponding three local rjs quasi-spin operators q l p. such that 
I = ±,z obey a SU{2) algebra and have the following ex- 
pression in terms of rotated-electron operators, 

qi = ) f ; % = (n?M - 1/2) • (55) 

Here = ti% ± i qi where x.y.z denotes the Cartesian 
coordinates. The relation of these r\s quasi-spin operators 
to the original electron creation and annihilation opera- 
tors involves the unitary transformation of Eq. (|43|) . 

Within the present rotated-electron operational formu- 
lation, three related elementary objects naturally emerge 
that make the model's global symmetry explicit. The 
operators provided in Eq. (|50[) create and annihilate 
spin-less and 77-spin-less c fermions whose local density 
operator, Eq. (|5l"j) , is directly related to the generator 
of the global c hidden i7(l) symmetry, as given in Eq. 
(|5"3"]l . The c fermions carry the charges of the rotated 
electrons that singly occupy sites. Moreover, the three 
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rotated local spin operators si s and the three rotated 
local ?7-spin operators si such that I = ±, z given in 
Eq. (|54p are associated with the spin-1/2 spinons and 
7?-spin-l/2 77-spinons, respectively, as defined here. The 
spin-1/2 spinons carry the spin of the rotated electrons 
that singly occupy sites. The c fermion holes describe 
the degrees of freedom associated with the c hidden U(l) 
symmetry of the sites doubly occupied and unoccupied 
by the rotated electrons. The 77-spin degrees of freedom 
of these sites are described by the 77-spin projection —1/2 
77-spinons (rotated-electron doubly occupied sites) and 77- 
spin projection +1/2 ?7-spinons (rotated-electron unoccu- 
pied sites). 

Within our representation, the local operators s?. c , 
sh, c , and si a where / = ±, z and a = s,rj can be ex- 
pressed in terms of only the c fermion local density op- 
erator g£ given in Eq. (|5 1 [) and three local rjs quasi-spin 

operators qL of Eq. ([55} as follows, 

S f ) j,C = Qrj J ^fj,C = 0- ~ I 

I = ±,z. (56) 

The expressions of the local spinon operators si g and 
local 77-spinon operators si ^ provided here are a con- 
firmation that the corresponding spin SU(2) and 77-spin 
SU (2) symmetries are particular cases of the rjs quasi- 
spin SU(2) symmetry. Specifically, they are associated 
with the SU (2) algebra representations involving the (i) 
spin-up and spin-down rotated-electron singly occupied 
sites and (ii) rotated-electron doubly-occupied and un- 
occupied sites, respectively. Indeed, the c fermion and 
c fermion hole local density operators q'i and (1 — ) 
play in the expressions of these operators provided in Eq. 
(|56f the role of projectors onto such two sets of lattice- 
site rotated-electron occupancies, respectively. 

The relations given in Eq. ([56)) for the operators si s 

and si 7] are equivalent to the following expression of the 
local 77s quasi-spin operators gl in terms of those of the 
former operators provided in Eq. (|54p . 

4 =4,,« + 4^. l = ±,z- (57) 

We emphasize that the c fermion operators, Eq. (|5(jp. 
and the spinon and 77-spinon operators defined by Eqs. 
(|54|) . (f55|) . and |56|) are mapped from the rotated-electron 
operators by an exact local unitary transformation that 
does not introduce constraints. Given their direct re- 
lation to the generators of the model's extended global 
symmetry, their occupancy configurations naturally gen- 
erate representations of the corresponding global sym- 
metry algebra. Consistent with the lack of constraints 
of such a local unitary transformation, inversion of the 
relations given in Eqs. ([5CT)l and (|55|) fully defines the 
rotated-electron operators in terms of the c fermion and 



77s quasi-spin operators as follows, 

«*,t = M,c - + /t >c Q + ql ) , 

%4 = (/^c + e^'/^o)^- (58) 

As given in Eq. (|A11[) of Appendix [A] that the c fermion 
operators commute with the 77s quasi-spin operators is 
behind the form of the expressions given here, whose c 
fermion creation and annihilation operators are located 
on the left-hand side. 

The c fermion operator and 77s quasi-spin operator ex- 
pressions in terms of rotated-electron creation and anni- 
hilation operators given in Eqs. ([50]) and (f55|) . respec- 
tively, are except for unimportant phase factors similar 
to those considered in the studies of Refs^r— in terms 
of electron creation and annihilation operators. Our op- 
erational representation has the advantage of rotated- 
electron single and double occupancy being good quan- 
tum numbers for all finite interaction values. On the 
other hand, the operator expressions provided in Eqs. 
([50)1 and ([55)1 differ from those of Refsii^ by unimpor- 
tant phase factors. 

Since for finite U/t values the Hamiltonian H of Eq. 
(TTJ) does not commute with the unitary operator V = 

e~ s , when expressed in terms of the rotated-electron cre- 
ation and annihilation operators of Eq. (|43D it has an 
infinite number of terms, 

H = VHV f = H+[H,S] 

+ ±[[H, §],§] + .... (59) 

The commutator [H,S] does not vanish except for U/t—> 
00 so that H 7^ H for finite values of U/t. 

Provided that both U / 1 is finite and one accounts for all 
higher-order terms on the right- hand-side of Eq. (|59[) , the 
corresponding expression refers to the Hubbard model. 
This is in contrast to the physical problem studied in 
Refs.— ~— , for which the rotated creation and annihila- 
tion operators of Eq. ([4"3")l refer to electrons. Thus ex- 
cept for U/t — > 00 within the physical problem studied in 
Refs<2Ir— the Hamiltonian given in Eq. ([59")) is not the 
Hubbard Hamiltonian. Instead, it is a rotated Hamilto- 
nian for which electron double occupancy and single oc- 
cupancy are good quantum numbers. On the other hand, 
for the alternative physical problem studied here and in 
Refs.— ^ the rotated creation and annihilation operators 
of Eq. (|43[) refer to rotated electrons and the Hamilto- 
nian provided in Eq. (|59[) is the Hubbard Hamiltonian. 

The latter Hamiltonian may be developed into an ex- 
pansion whose terms can be written as products of the 
rotated kinetic operators T 7 given in Eq. (|A3[) of Ap- 
pendix [X] where 7 = 0, ±1. The corresponding order of 
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a given Hamiltonian term refers to the number of such 
rotated kinetic operators T 7 independently of their type, 
7 = 0, ±1. To fourth order such an Hamiltonian reads, 

H = A® + ff« + + ff(3) + fiW + ... , 
#(°) = UV D ; H {1) =tf , 



u 



T_iT. 



l-t+i j 



+3 i 

R{3) = 77a[T-iT f + i--(r_if + ifo + f f_if + i)] 

(7 2 

t 4 

i? (4) = [f-lfof+ifo + f f_if T+l 

+ f_iT + if_if+i - ^(f^f+ifo 2 + f 2 f-if +1 ) 
+ (2f f_if + if - T-xf+if 2 - f 2 f_ 1 f +1 )] , 



Here, 



— real — number parameter . 



(60) 



(n^ it -l/2) (n r -,;-l/2) 



(61) 



is the rotated-electron interaction operator. That it ap- 
pears only once in the Hamiltonian expansion whose 
leading-order terms are given in Eq. (|60p follows from 
the derivation of that expansion systematically using the 
commutator. 



[Vb,T 7 ] = 7IL 



7 = 0,±1 



(62) 



We recall that except for U/t — > oo one has that 
f ^ f , f+i f+i, and T_i ^ T-\. Expressing the 
Hamiltonina expression of Eq. in terms of electron 
creation and annihilation operators gives for large U/t 
values a similar expansion. However, for the intermedi- 
ate U/t values of interest for our study the few first terms 
of the Hamiltonian t/U expansion given in of Eq. (|60|) in 
terms of rotated-electron operators contain much more 
complicated higher-order terms when expressed in terms 
of electron creation and annihilation operators. 

That Hamiltonian expansion may be expressed in 
terms of the c fermion and rjs quasi-spin operators. This 
is achieved by combining the rotated-electron operator 
expressions provided in Eq. with those of the 

rotated-electron interaction operator given in Eq. (|61j) 
and three rotated kinetic operators Tq, T_i, and T + \ pro- 
vided in Eq. (|A"3f of Appendix [S] 

If a rotated-electron term of an operator expansion in 
terms of rotated-electron creation and annihilation oper- 
ators does not preserve the numbers of rotated-electron 
singly and doubly occupied sites, we call it off-diagonal. 
An interesting technical detail is that up to third order all 
diagonal terms of the Hamiltonian expression provided in 
Eq. (j6"0")) are generated by the leading-order term of the 



operator S, which is given in Eq. (|4"4"| . Indeed, when 
expressed in terms of electron operators the Hubbard 
Hamiltonian provided in Eq. (JTJ) does not contain any 
off-diagonal terms with more than two electron opera- 
tors. (In this case the off-diagonal terms are electron off- 
diagonal terms, which refer to electron doubly occupied 
sites.) 

Only the Hamiltonian terms ffW, and 

iT 3 * 1 to third order given in Eq. (|6"0")) are universal. In- 
deed, the form of the terms of fourth and larger order 
is different for each electron - rotated-electron unitary 
transformation. For the fourth-order term iT 4 ) given in 
that equation only the real-number parameter 9 value is 
not universal, being unitary-transformation dependent^ 2 -. 
For instance, the methods of Refs^ and^i refer to two 
different electron - rotated-electron unitary transforma- 
tions whose 9 values are 9 = and 9 = 1/4, respec- 
tively. Moreover, one of the methods of Refi^ 2 - refers 
to an electron - rotated-electron unitary transformation 
whose 9 value is 6 = 1/2. Its value for the electron 
- rotated-electron unitary transformation whose unitary 
operator is uniquely defined by the matrix elements of 
Eqs. (|46|) and (|47|) remains an open issue. Fortunately, 
these Hamiltonian terms multiplying the parameter 9 
vanish at half filling so that this does not affect the en- 
suing section studies. 

The non-universal Hamiltonian terms are all reducible 
with respect to the subspaces with fixed values of rotated- 
electron single and double occupancies. That is, they 
contain hopping processes that do not originate from ex- 
citation between these subspaces, T)T_iT + iTo, nor ter- 
minate once a rotated-electron or rotated-hole is returned 
to a subspace with larger singly-occupancy, for example, 
T_iT + iT_iT + i. All these processes can be viewed as 
arising from the specific transformation V r ^|4'; ri ; ?)si00 ) of 
the U/t — > oo energy eigenstates within the subspaces 
with fixed values of rotated-electron single and double 
occupancies. Thus the infinite electron - rotated-electron 
unitary transformations differ in the processes within 
each subspace with fixed values of these occupancies. 



VI. GENERAL S 3 = 1 SPIN SPECTRUM 
WITHIN THE SPINON REPRESENTATION 

As discussed in Sec. U the usual spin- wave theory 
does not describe the neutron scattering of LCO. Here we 
study the S s = 1 spin-triplet spectrum of the half-filled 
Hubbard model on the square lattice by means of the 
spinon representation that emerges from the above more 
general c fermion and r/s quasi-spin operator formulation, 
which is that suitable for the LCO intermediate interac- 
tion range U/t £ (6,8). (In units of the bandwidth, 8t, 
this gives U/8te (0.75,1).) 

For very large U/t values the Hubbard model may be 
mapped onto a spin-only problem whose spins are those 
of the electrons that singly occupy sites. However, for 
intermediate U/t values electron single occupancy is not 
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a good quantum number so that such a mapping breaks 
down. On the other hand, the rotated electrons of our 
operator representation have been constructed to make 
rotated-electron single and double occupancy good quan- 
tum numbers for U/t > 0. This is why our spinons are 
well defined for the LCO intermediate interaction range 
U/t E (6,8). Indeed they are the spins of the rotated 
electrons that singly occupy sites. In the laige-U/t limit 
the rotated-electrons become electrons, so that one re- 
covers the known standard results. 

Within our operator formulation, the Hubbard model 
in the vanishing rotated-electron doubly occupied sites 
number and unoccupied sites number subspace (VDU 
subspace) can be mapped onto a spin-only problem for 
all U jt finite values. In the VDU subspace the num- 
ber of spin-1/2 spinons equals that of rotated-electrons, 
electrons, and sites N = N a . Since there are no rotated- 
electron doubly occupied or unoccupied sites there are no 
?y-spinons. Hence the number of //-spin SU (2) symme- 
try algebra representations vanishes and that symmetry 
does not play any role. Furthermore, although there are 
TV = N a c fermions, their c momentum band associated 
with the operators of Eq. ([52")) is full. Hence the degrees 
of freedom associated with the c fermion occupancy con- 
figurations that generate the c hidden U(l) symmetry al- 
gebra representations are frozen and the Hubbard model 
in the subspace under consideration may be mapped onto 
a spin-only problem, as confirmed below. 

For U/t — > oo the VDU subspace is the only one for fi- 
nite excitation energy. For the finite-?/ /t spin excitations 
that preserve the electron number N = N a considered 
in the following, it is the only subspace within a finite 
excitation-energy window, uj € (0, 2Amh)- Here 2Amh 
is the Mott-Hubbard gap. Below we calculate its U/t 
dependence for the LCO intermediate interaction range 
U/t e (6,8) by DMRG. Our goal is to check whether 
the relevant spin energy spectrum that emerges from our 
VDU subspace spin-only problem is indeed contained in 
the excitation-energy domain u € (0, 2Amh) f° r which 
it is valid. 
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Figure 3: The Mott-Hubbard gap 2A M h DMRG numerical 
results on two different width cylinders along with the ap- 
proximate analytical expression of Eq. (|63[) (solid line) as 
a function of U/t. The DMRG points seem to be consis- 
tent with for the half-filled Hubbard on the square lattice the 
Mott-Hubbard gap being finite for U/t > and vanishing in 
the U/t -> limit. 



Our DMRG calculations refer to the single-particle 
gap. They have been performed both for 10 x 4 and 
10 x 6 Hubbard cylinders. The chemical potential was 
set to U/2 and two states were targeted, one with N 
particles and the other with N — 1. (Targetting N + 1 
electrons would have given the same results.) Both states 
were put into the same density matrix in the traditional 
multi-state targeting DMRG approach. Thus, the same 
truncation error applied to both states, leading to signif- 
icant error cancellation. The resulting gap at each sweep 
was plotted versus the maximum truncation error in the 
sweep, yielding approximately linear behavior, and allow- 
ing the extrapolation to zero truncation error. The error 
estimate is roughly the size of the extrapolation from the 
last point. From 1800 (10 x 4 Hubbard cylinder) to 6000 
(10 x 6 Hubbard cylinder) states were kept. 

Here we report the corresponding magnitudes of the 
Mott-Hubbard gap 2A MH . For the range U/t e (4,20) 
we find that, 



A. The energy range of our spin-only quantum 
problem 

From the interplay of the model's symmetries with our 
operator formulation that makes these symmetries ex- 
plicit, one straightforwardly confirms that the minimum 
energy for creation of one rotated-electron doubly occu- 
pied site or one rotated-electron unoccupied site at fixed 
electron number N = N a onto the n = 1 and m = 
ground state is indeed given by the Mott-Hubbard gap, 
2Amh- Its magnitude is twice that of the single-particle 
gap, Amh- In order to define the energy range of our 
study, here we calculate the Mott-Hubbard gap 2Amh 
dependence on U ft for a domain containing the LCO 
range U/t £ (6,8). 



2A 



AIH 



U 



tanh 



U/lt 



Ve-y+y/Uf-yt 



6 



6 ' 
(63) 

gives quantitative agreement with our numerical DMRG 
calculations for the Mott-Hubbard gap 2Amh depen- 
dence on U/t. The DMRG points for that gap are plotted 
in Fig. [3] along with the curve obtained from the approx- 
imate analytical expression, Eq. ()63[) . 

For instance, our DMRG calculations for 10 x 6 Hub- 
bard cylinders give 2A M h ~ 2.78(4) t for U/t = 6 and 
A M h = 4.30(4)4 for U/t = 8. This leads to a range 
2 Amh G (816 meV, 1442 meV) for U/t G (6, 8). Here we 
used the t magnitudes t ps 293 meV and t ps 335 meV 
for which the model describes the LCO neutron scat- 
tering for U/t = 6 and U/t = 8, respectively. For the 
U/t = 6.1 value used in some of our calculations, we 
find A M h ~ 2.81(0) t from the DMRG analysis, so that 
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A M h ~ 829 meV for t « 295 meV. 

Optical experiments overestimate the charge-transfer 
gap magnitudes of the parent insulating compounds^. 
On the other hand, by measuring the Hall coefficient Rh 
in LCO, the studies of Ref^ have estimated the energy 
gap over which the electron and hole carriers are ther- 
mally activated, which corresponds to the Mott-Hubbard 
gap, to be 2Amh ~ 890 meV. Remarkably, this magni- 
tude is within the range 2/S.mh £ (816 meV, 1442 meV) 
of our above theoretical predictions for U/t € (6,8). 
Our theoretical approach based on the combination of 
our DMRG results with the U and t values for which 
agreement with the LCO neutron-scattering agreement 
is reached leads to 2Amh ~ 890 meV for U/t w 6.3. 
Below we consistently confirm that the spin-triplet ex- 
citation spectrum calculated for the Hubbard model in 
the VDU subspace is contained in the energy window 
ui e (0, 2Ajv/#) found here. 

Note that the DMRG points of Fig. [3] seem to be 
consistent with the Mott-Hubbard gap being finite for 
U/t > and vanishing in the U/t — > limit. 

B. The Hubbard model in the VDU subspace 

Let us confirm that within our operator representation 
the half-filled Hubbard model on the square lattice in 
the VDU subspace can for U/t > be expressed solely 
in terms of spinon operators. Indeed, accounting for the 
lack of both rotated-electron doubly occupied sites and 
unoccupied sites, upon writing the Hamiltonian of Eq. 
([60| in the VDU subspace, one finds that all its terms 
of odd order vanish and the terms of even order given in 
that equation simplify to, 

= UV C , 

H (2) = ~f_if + i, 

t 4 

Hi4) = jp[T-iT+iT-iT +1 

l f 3^ _ f_ lf 2f +l] _ (64) 

We have then expressed the Hamiltonian terms of even 
order as those provided in Eq. (|64l) in terms of the c 
fermion and rjs quasi-spin operators. This has been done 
by combining the rotated-electron operator expressions 
provided in Eq. ([55)1 with those of the three rotated ki- 
netic operators To, T!_i, and T+i given in Eq. (|A3|) of 
Appendix |XJ Since the states that span the VDU sub- 
space are generated only by rotated-electron singly occu- 
pancy configurations, the projectors g£. and (1 — gS.) in 
the expressions of Eq. (|5t?|) can be replaced by the cor- 
responding eigenvalues 1 and 0, respectively. One then 
finds that s^. = q 1 -. in the VDU subspace, so that the 
77-spinon operators do not play any role. Hence in it the 
77s quasi-spin operators q^, reduce to the corresponding 

spinon operators sL „, where I = ±, z. 



Moreover, after some algebra involving the anti- 
commutation and commutation relations given in Eqs. 
(|ATB - ([A13| of Appendix E] one finds that all contribu- 
tions involving the c fermion creation and annihilation 
operators can be expressed only in terms of local oper- 
ators <jj. In the VDU subspace one can then replace 
these operators by their eigenvalue 1. Thus the Hamilto- 
nian terms of Eq. (|64[) can be expressed only in terms of 
spinon operators. Importantly, this holds as well for the 
remaining Hamiltonian terms of higher even order omit- 
ted in that equation. Moreover, all Hamiltonian terms 
of odd order vanish and the zeroth-order term becomes 
a mere constant, H(") = [U/i] N a , and may be ignored. 
The Hamiltinonian terms of second and fourth order of 
Eq. (|64|) may after some algebra then be rewritten as, 



H [2) = Jj E ^s-^-l], (65) 

(hh) 

and 




(^fjjjS ' Srj 3 ,s)(Srj 2 ,s ' ^r 3 ' 4 ,s)j > (66) 



respectively. Here the spinon operator s? jtS has opera- 
tor Cartesian components s5 . = hist + sZ 1, s v - = 

*r rj.s 2 L rj,s 1 rj.s-l' rj,s 

jr[si. s — s~. s ], and st s and refers to the spin of a ro- 
tated electron that singly occupies the site of real-space 
coordinate fj. The spinon operators s^. s and si s are 
those given in Eq. (|54[) . Furthermore, in the expres- 
sions of Eqs. (|65p and (|66p the summation (j'ij'2) runs 
over nearest-neighboring sites and Djji = 1 for the real- 
space coordinates fj and fj' corresponding to nearest- 
neigboring sites and Djji = otherwise. 

For very large U/t values when electron single and dou- 
ble occupancy become good quantum numbers and thus 
the rotated electrons become electrons the spinon opera- 
tors Sf. <8 become the usual spin operators Sf jjS and Eqs. 
(|65|) and (|66|) recover the corresponding spin-only Hamil- 
tonian terms obtained previously by other authors^. On 
the other hand, for the intermediate U/t values of inter- 
est for LCO the terms of the Hamiltonian t/U expansion 
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given in of Eq. in terms of spinon (rotated-electron) 
operators contain much more complicated higher-order 
terms when expressed in terms of electron creation and 
annihilation operators. 

C. The absolute ground state of the Hubbard 
model on the square lattice 

The antiferromagnetic long-range order of the half- 
filled Hubbard model on the square lattice ground state 
follows from a spontaneous symmetry breaking mecha- 
nism that occurs in the thermodynamic limit N a — > oo. 
It involves a whole tower of low-lying energy eigenstates 
of the finite system. They collapse in that limit onto the 
ground state. 

Importantly, both that ground state and the excited 
energy eigenstates that collapse onto it as iV a — > oo be- 
long to the VDU subspace. One may investigate which 
energy eigenstates couple to the exact finite N a >• 1 and 
n = 1 and m = ground state \GS) via the operator, 

^^FB^'i" l = ±,z. (67) 
' Vq 3=1 

We insert the complete set of energy eigenstates as fol- 
lows, 

(GS\(m l s ) 2 \GS) = (GS\m l s \y lrJvs , u ) 

= J2 \( G S\™ l s\*ir^,u)\ 2 ! l = ±,X3. (68) 

Only energy eigenstates \^i r ,i with excitation mo- 
mentum k = tt and quantum numbers S v = 0, 2S C = 
N a = N, S s = 1, and S* = 0, ±1 contribute to the sum 
of Eq. (|68|) . We recall that the quantum numbers S 71 = 
and 2S C = N a = N remain unchanged and thus are the 
same as those of the ground state \GS). We denote by 
|*it) the S s = 1,S V = 0, 2S C = N a = N, and k = [tt, tt] 
lowest spin-triplet state whose excitation energy behaves 
as 1/N a for finite N a > 1. For the range U/t > 4 of 
interest for our studies the contribution from this lowest 
spin triplet state is by far the largest. For instance, for 
the related spin- 1/2 Heisenberg model on the square lat- 
tice the matrix-element square |(GS'|m / s | v I'iT)| 2 exhausts 
the sum in Eq. (|6"8"]l by more than 98.7%2I. A similar be- 
havior is expected for the Hubbard model on the square 
lattice, at least provided that U/t > 4. 

The special properties with respect to the lattice sym- 
metry group of the lowest energy eigenstates contribut- 
ing to the linear Goldstone modes of the corresponding 
S s = 1 spin-wave spectrum reveal the space-symmetry 
breaking of the N a — > oo ground state. In the present 
case of the half-filled Hubbard model on the square lat- 
tice the translation symmetry is broken. Hence as found 



here both the k = [0, 0] and k = [tt, tt] excitation mo- 
menta appear among the lowest energy eigenstates con- 
tributing to the linear Goldstone modes of the S s = 1 
spin- wave spectrum. However, that the transitions to the 
lowest spin-triplet state I'J'it) of momentum k = [tt,tv] 
nearly exhaust the sum in Eq. (|68|) is consistent with the 
first- moment sum rules of an isotropic antiferromagnet, 
such that no weight is generated by states of momentum 
k= [0,0]. 

One of the few exact theorems that apply to the half- 
filled Hubbard model on a bipartite lattice and thus on 
a square lattice is that for a finite number of lattice sites 
N a its ground state is a spin-singlet stated. The stud- 
ies of Refsjiii^ use an operator representation that dif- 
fers from that used here only by unimportant phase fac- 
tors. Such studies provide evidence that the n = 1 and 
m = ground state is the only model's ground state that 
is invariant under the electron - rotated-electron unitary 
transformation. For N a ^> 1 the results of those refer- 
ences reveal that its spin-singlet configurations refer to 
N a /2 = N/2 independent spin-singlet two-spinon pairs. 
Most of the weight of these spin-singlet two-spinon pairs 
stems from spinons at nearest-neighboring sites yet they 
have finite contributions as well from spinons located at 
larger distances. 

Our spinon representation has been constructed to 
make such N/2 spin-singlet spinon pairs correspond to 
spin-neutral objects that obey a hard-core bosonic alge- 
bra. One can then perform an extended Jordan- Wigner 
transformation that maps them onto N/2 si fermionsH. 
(In the index si the number 1 refers to one spin-singlet 
spinon pair.) The corresponding si fermion momentum 
band is full for the n = 1 and m = absolute ground 
state. It has a momentum area 2tt 2 and coincides with 
an antiferromagnetic RBZ whose momentum q compo- 
nents obey the inequality, 

\q x \ + \q y \<lT. (69) 

As a result of its invariance under the electron - 
rotated-electron unitary transformation, the n = 1 and 
m = absolute ground state is the only ground state that 
for U/t > belongs to a single V tower. Hence both for 
it and its spin-triplet excited states that belong to the 
VDU subspace the si boundary-line momenta qssi are 
independent of U/At. Consistent with Eq. their 
Cartesian components qssix and qBsiy obey the equa- 
tions, 

iBslx ± qBsly = TT , 

or q Bs \x ± QBsly = -tt ■ (70) 

Hence the si boundary line refers to the lines connecting 
[±tt,0] and [0,±tt]. 
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D. The spin excitations and the ground-state 
spinon d-wave pairing 

Within our spinon operator representation the S s = 1 
spin-triplet excitations relative to the n — 1 and m = 
absolute ground state involve creation of two holes in the 
si band along with a shift Tv/N a of all discrete momentum 
values of the full c band. Under such an excitation one 
of the N a /2 — N/2 spin-singlet spinon pairs is broken. 
This gives rise to two unbound spinons in the excited 
state whose three occupancy configurations generate the 
three spin-triplet states of spin projection S* = 0, ±1. 
In the case of such spin-triplet excitations the occupancy 
configurations of the two holes arising in the si fermion 
momentum band may simulate the motion of the two 
unbound spinons relative to a background of N/2 — 1 
spinon pairs, or vice versa. 

The general spin-triplet spectrum has within the 
present spinon representation the following form, 



oj(k) 
k 



-e s i(g) - 
'-q-q 



(<?') 



(71) 



where if = ±[7r, ±7r], q and q' are the momentum values 
of the emerging two si fermion holes, and e s i{q) is the 
corresponding si fermion energy dispersion. Indeed, the 
results of RefJ^ provide evidence that for the Hubbard 
model on the square lattice in the VDU subspace the 
si fermion momentum q is a good quantum number, so 
that one can define a corresponding energy dispersion. 
However, in contrast to ID this property does not hold 
for the more general problem of that model in its full 
Hilbert spaced. 

The Hubbard model on the also bipartite ID lattice has 
the same extended global symmetry than on the square 
lattice. Hence for it an operator representation similar 
to that used here may be introduced. The exact Bethe- 
anstaz solution then implicitly performs the summation 
of all Hamiltonian terms of even order whose leading- 
order terms are given in Eqs. (|65p and (|6"6"|) . This leads 
to a si fermion band e s \(q) that in the U/t — > limit 
equals the occupied part of the electron non-interacting 
dispersio n 38 i 39 . The main effect of increasing the U/t 
value is decreasing the si fermion band e s i(q) energy 
bandwidth. It decreases from 2t as U/t — > to zero for 
U/t -> oo. 

As discussed above, expression of the Hamiltonian in 
terms of rotated-electron operators leads for the interme- 
diate U/t € (6, 8) range to a quantum problem in terms 
of rotated-electron processes similar to the correspond- 
ing large-C/ /t quantum problem in terms of electron pro- 
cesses. At half filling the main effect of decreasing U/t 
is the increase of the energy bandwidth of an effective 
band associated with the spinon occupancy configura- 
tions. Such an effective band is the si energy dispersion. 
Consistent with and partially motivated by the exact ID 
results yet accounting both for the corresponding com- 
mon global symmetry and different physics, the rotated- 
electron studies of Refsii 3 - provide evidence that for the 




Figure 4: The energy-momentum space limits of the spin S 3 = 
1 excited states spectrum of Eq. (|71[) for U/t = 6.1, t = 
295 meV, and k x and k y in units of 2n. States whose energy 
is for a given k lower than that of the intermediate spin-wave 
sheet as well as those of any energy and equivalent momenta 
[0, 0] = [0, 2tt] = [2tt, 0] = [2tt, 2tt] do not contribute to the 
spin spectral weight. 



model on the square lattice the effective si energy dis- 
persion e s ±(q) involves an auxiliary dispersion, 

W° 

4 (0) = f- [cos q x + cos q y ] . (72) 

In the U/t — > limit such an auxiliary dispersion reaches 
its maximum energy bandwidth. Similarly to ID, in that 
limit it is expected to become the occupied part of the 
electron non-interacting dispersion. The main effect of 
increasing U/t is to decrease the energy bandwidth of 
that dispersion and thus the magnitude of the energy 
scale Wgi in Eq. (|72")l . so that for half filling it changes 
from = At as U/t -> to = for U/t -+ oo. 

However, the n = 1 and m = ground state of the ID 
half-filled Hubbard model has no antiferromagnetic long 
range order as N a — > oo. In the presence of that order, 
provided that U/t is not too small so that one can ignore 
the amplitude fluctuations of the corresponding order pa- 
rameter, the problem can be handled for the model on 
the square lattice by a suitable mean-field theory. Within 
it the occurrence of that order is described by a si energy 
dispersion of the general formii 3 -, 



esi(q) = -J\e° al (q)\ 2 + \^si(q)\ r ' 



(73) 



Here £®i{q) is the auxiliary dispersion given in Eq. (|73p 
and the gap function |A s i(g)| is to be determined from 
comparison with the spin-triplet spectrum obtained from 
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the standard formalism of many-body physics by sum- 
ming up an infinite number of ladder diagrams. (We 
note that as explicitly shown in Ref£ the RPA studies of 
Sec. HVI are equivalent to summing up an infinite number 
of such diagrams.) 

We profit from symmetry and limit our analysis of the 
spin spectrum of Eq. (|7T|) to the sector k x G (0,7r) and 
k y G (0, k x ) of the (k,u>) space. Surprisingly, quantitative 
agreement with the results obtained from summing up an 
infinite number of diagrams is reached provided that the 
si dispersion gap function refers to a ci-wave si fermion 
spin-singlet spinon pairing, 



|A s i(?)| 



2 



cos q x - cos q y | 



(74) 



Moreover, from comparison with many-body physics re- 
sults one finds that the inelastic coherent spin- wave spec- 
trum is generated by processes where q points in the 
nodal direction and q' belongs to the boundary of the 
si band reduced zone. The remaining choices of q and q' 
either generate the inelastic incoherent continuum spec- 
tral weight or vanishing weight, respectively. 

For this choice of the momenta of the two emerging si 
fermion holes one finds from the use of Eqs. (|72l) - ([74|) 
that the spin-wave spectrum corresponds to a surface of 
energy and momentum given by, 



«® = y 



k x + kx 



sin 



k x kr 



q-q 



(75) 



This is a particular case of the general spin spectrum 
of Eq. (|7ip , which refers to the following choices of the 
momenta ff, q, and q 1 , 
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for the sub-sector such that k x G (0,7r), k y G (0, k x ) for 
k x < t/2, and k y G (0, tt — k x ) for k x > ir/2. Moreover, 
for the sub-sector such that k y G (0,7r), k x G (tt — k yi ir) 
for k y < 7r/2, and k x G (k y ,Tt) for k y > 7r/2, respectively, 
it corresponds to the following choices of the momenta tt, 
q, and q' , 
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Note that the components of the si band momenta q 
appearing in Eqs. (|76[) and (|77|) are such that q x — q y = 
—tt and thus belong to the half-filling si boundary line 



defined by Eq. ([TT?)) , whereas those of the momenta q' in 
the same equations obey the relation q' x 
point ir x1 J ~' J ' — i: 



q' so that 
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Figure 5: The theoretical spin spectra, Eqs. (|B1[) - (|B9[) . (solid 
lines) plotted in the second BZ for U/t » 6.1 and t ~ 0.295 
eV and thus f/P — 565.6 meV and W°\ = 49.6 meV and the 
experimental data of Ref.— (circles) in meV. The momentum 
is given in units of 2tt. From Ref.—. 

For the values U/t = 6.1 and t = 295 meV used in 
Sec. IIVI in our study of the LCO spin spectrum one finds 
that « t/5.95 and fi° w t/0.5216 in the expressions 
of Eqs. ((72J), (O, and f^J), so that w 49.6 meV 
and ^° f=a 565.6 meV. The spin- wave spectrum of Eq. 
([75)1 calculated for these and /i° values refers to the 
middle surface plotted in Fig. |4j Its expressions corre- 
sponding to the high symmetry directions in the BZ are 
given in Appendix [Bj The corresponding curves are plot- 
ted in the top panel of Fig. [21 along with those obtained 
from the many-body physics by summing up an infinite 
number of ladder diagrams and the LCO experimental 
points of Ref»2. In Fig. [5] the same curves arc plotted 
together with the LCO experimental points of RefA We 
emphasize that the intermediate sheet plotted in Fig 01 
which corresponds to the general spin-wave spectrum of 
Eq. ([75)1 , also fully agrees with the results of the experi- 
mental studies reported in Refsi^. 

The studies of Ref>£ are limited to the spin- wave spec- 
trum. Following the agreement of the spin-wave spec- 
trum of Eq. (|75p obtained from the general spin spec- 
trum of Eq. (|7ip with both results from the many-body 
physics and LCO neutron-scattering experimental points 
of Refs<2i^ here we consider it for all choices of the si 
fermion hole momenta q and q'. The corresponding 
energy-momentum space domain of the spin S s = 1 ex- 
cited states whose spectrum is provided in Eq. ([7T|) is 
represented in Fig. 0] for U/t = 6.1 and t — 295 meV. 
A similar spectrum is obtained for the values U/t = 8.0 
and t = 335 meV of Ref£. 

The largest energy of the general spin spectrum of the 
half-filled Hubbard model on the square lattice in the 
VDU subspace represented in Fig. @] is 566 meV. Hence 



the whole spin spectrum of Eq. (|71|) represented in that 
figure is contained in the energy window cj G (0, 2Amh) 
of the corresponding VDU subspace. Indeed, from com- 
bination of the results of our DMRG calculations with 
the t magnitudes that lead to agreement with the spin- 
wave spectrum of LCO we have found that 2Amh G (816 
meV, 1442 meV) for U/t G (6, 8). 

As mentioned above, the intermediate sheet of the 
general spin spectrum represented in Fig. U refers to 
the spin-wave spectrum. For each excitation momen- 
tum k, states of energy lower than the latter spectrum 
do not contribute to the form- factor weight. Further- 
more and consistent with the first-moment sum rules of 
an isotropic antiferromagnet, no and nearly no weight 
is generated by states of any energy and momentum 
[0,0] = [0,2tt] = [2tt,0] = [27r,27r] and near it, respec- 
tively. 

Unfortunately, in its present form our spinon-operator 
method does not provide the detailed continuum weight 
intensity distribution. However, it is expected that, simi- 
larly to the Heisenberg model case^, its energy-integrated 
intensity follows the same trend as the spin-wave inten- 
sity. Analysis of Fig. 0] reveals that for momentum [n,0] 
there are no excited states of energy higher than the spin 
waves. Thus at momentum [7r,0], the continuum weight 
distribution energy-integrated intensity vanishes or is ex- 
tremely small, due to si band four-hole processes. Given 
the expected common trend of both intensities, this is 
consistent with a damping of the spin-wave intensity at 
momentum [7r,0], as observed in the recent high-energy 
inelastic neutron scattering experiments of Ref.— but not 
captured by the Fig. [5] RPA intensity. Since we find 
that for the half-filled Hubbard model on the square lat- 
tice the continuum weight distribution energy-integrated 
intensity vanishes or is extremely small at momentum 
[7r,0], we predict that for it the corresponding spin- wave 
intensity at momentum [n, 0] is also damped. Hence we 
expect that such a behavior is absent in Fig. [5]due to the 
RPA used in the calculation of the spin-wave intensity. 
However, for all other momentum values the RPA results 
are expected to be a quite good estimate of the model's 
spin-wave intensity. 

The d-wave spin-singlet spinon pairing that follows 
from the energy gap Eq. (|74[) emerged here from im- 
posing quantitative agreement with the spin-wave spec- 
trum obtained from summing up an infinite number of 
diagrams. We emphasize that such a type of spinon 
pairing is not inconsistent with the ground-state anti- 
ferromagnetic order provided that the weights of the cor- 
responding spinon pairs fall off as a power law of the 
spinon distance whose negative exponent has absolute 
value smaller than 5. This was confirmed in Ref4£ for 
spins associated with electrons yet holds as well for the 
present spinons, which refer to the spins of the rotated- 
electrons that singly occupy sites. We emphasize that 
such a pairing does not refer to electrons or rotated elec- 
trons. For the n = 1 and m = absolute ground state it 
corresponds to N/2 spinon pairs, which describe the spin 
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U/t 


6.1 


6.5 


8.0 


10.0 


W T 


1.643 


1.671 


1.762 (1.7782) 


1.848 (1.8462) 


Wsw 


0.808 


0.799 


0.761 


0.714 


W<450 


1.571 


1.593 


1.663 


1.730 


W / >450 


0.072 


0.078 


0.099 


0.118 



Table II: Several spectral weights in units of fi% as denned in 
the text for several U/t values and some results from Refs.— i— . 

degrees of freedom of the N rotated electrons. Indeed, 
for that ground state all rotated electrons singly occupy 
sites. 



VII. COMPARISON OF THE PREDICTED 
SPECTRAL WEIGHTS WITH THOSE IN THE 
LCO HIGH-ENERGY NEUTRON SCATTERING 

As discussed in Sec. IIV1 the total spin-weight sum-rule, 
fi B 2(1 — 2d), of Eq. Q can in units of \j? B be written as, 

W T = 2(1 - 2d) « Wsw + Wco + 4(m^f ) 2 . (78) 

Here Wsw = Z d [2(1 - 2d) - A{m A F) 2 ] is the integrated 
spectral weight associated with the spin-wave intensity, 
Eqs. (JMD and W CO = 0—Z d ) [2(l~2d)~4(m AF ) 2 } 
that of the remaining inelastic spin spectral weight as- 
sociated with the continuum distribution, and 4(m^J?) 2 
refers to the Bragg-peak elastic part of the spin spectral 
weight. 

In Table |H] we provide the results of our calculations 
for several integrated spin spectral weights (in units of 
H 2 B ). This includes the total spin spectral weight Wt and 
the spin- wave intensity coherent spectral weight Wsw- 
In addition, in the table we provide the estimated mag- 
nitudes of the total spin spectral weight for excitation 
energy Huj <450 meV, IU< 4 5o = WW/0.71 + 4(ragf ) 2 , 
and the total spin spectral weight for excitation energy 
hu > 450 meV, IU>4 5 o = [Wt — IU<4 50 ]. The magni- 
tude of the spin spectral weight for excitation energy 
Hu> <450 meV, IU<45o, is derived by combining our theo- 
retical expressions with the observations of Ref^ that for 
the energy range up to about 450 meV, 71% and 29% of 
the weight corresponding to the inelastic response comes 
from the coherent spin-wave weight and incoherent con- 
tinuum weight, respectively. 

Our prediction for the magnitude of the spin spectral 
weight for excitation energy Huj <450 meV, TU<45o, varies 
between 1.6 for U/t « 6.1 and 1.7 (J? B for U/t « 8.0, 
in agreement with the experimental value 1.9 ± 0.3 fi B 
reported in Refj2. 

From our above analysis, the amount of spin spec- 
tral weight for excitation energy hu > 450 meV is small, 
H / >45o ~ 0.1 /j, b . By combining our spectral- weight re- 
sults with the boundaries of the spin-triplet spectrum 
plotted in Fig. [H such a small spin spectral weight is 
expected to extend to about 566 meV, mostly at and 
around the momentum [7T,7rj. 
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VIII. CONCLUDING REMARKS 

In this paper we have studied by means of the half- 
filled Hubbard model on the square lattice several open 
issues raised by the recent LCO neutron scattering data 
reported in Refi^. Our studies combined standard meth- 
ods such as RPA techniques involving a broken symmetry 
ground state with DMRG calculations on Hubbard cylin- 
ders and a spinon representation suitable to the LCO 
intermediate interaction range U/t £ (6,8). The latter 
emerges from a rotated-electron operator representation 
that has been constructed to ensure that rotated-electron 
single and double occupancy are good quantum numbers 
for U/t > 0. This assures that our spinons are well de- 
fined for the LCO intermediate U jt range. Indeed these 
spinons are the spins of the rotated electrons that singly 
occupy sites. In the large-U/t limit the rotated-electrons 
become electrons, so that one recovers the usual picture. 

Within this operator formulation, the Hubbard model 
in the VDU subspace considered here can be mapped 
onto a spin-only problem for the U /t range of interest for 
our studies. The spin excitations preserve the electron 
number. At fixed electron number the VDU subspace 
is the only one within a finite excitation-energy window, 
lu £ (0, 2A M h)- Here 2A M h is the Mott-Hubbard gap, 
whose U/t dependence we have calculated by DMRG. 
We have found that 2A MH € (816 meV, 1442 meV) for 
U/t € (6,8). Correspondingly, the largest energy of the 
general spin spectrum of the half-filled Hubbard model 
on the square lattice in the VDU subspace represented 
in Fig. [l]is 566 meV, so that it is fully contained in this 
energy window. 

The coherent part of the spin spectrum, which corre- 
sponds to the spin-wave spectrum, was complement ar- 
ily studied by a RPA analysis involving a broken sym- 
metry ground state and the spinon operator represen- 
tation. The former method was also used to calculate 
the spin-wave intensity momentum distribution. Quanti- 
tative agreement with both the spin-wave spectrum ob- 
tained by summing up an infinite number of ladder dia- 
grams and that observed in LCO neutron scattering ex- 
perimental studies is reached by the spinon method pro- 
vided that the initial n = 1 and m = ground state 
has d-wave spinon pairing. For that ground state such a 
pairing refers only to the rotated-electron spin degrees of 
freedom. Whether upon hole doping such a pre-formed 
d-wave spinon pairing could lead to rotated-electron d- 
wave pairing or even related electron d-wave pairing is 
an issue that deserves further investigations. 

Following the good quantitative agreement with the 
spin-wave spectrum, the spinon representation was used 
to derive the full S s = 1 spin-triplet spectrum represented 
in Fig. 2]for the U and t values suitable for LCO. From 
analysis of that figure we have found that for the mo- 
mentum [it, 0] there are no excited states of energy higher 
than the spin waves. Thus at momentum [it, 0], the con- 
tinuum weight distribution energy-integrated intensity 
vanishes or is extremely small. Such an intensity is ex- 



pected to follow the same trend as that of the spin waves. 
Hence this behavior is consistent with a corresponding 
damping of the spin-wave intensity at [it, 0] observed in 
the recent high-energy inelastic neutron-scattering exper- 
iments of Refi£. 

On the other hand, a resonant-inelastic x-ray scatter- 
ing study of insulating and doped La2- 2; Sr 2 ;Cu04 found 
a mode at 500 meV, at a momentum transfer [tt, 0)^. 
This 500 meV mode is observed only when the incident 
x-ray polarization is normal to the CuO planes. It could 
be a d-d crystal-field excitatio n 42 ' 43 , rather than a spin 
excitation. In case it is a spin excitation, one possible 
explanation given in Ref4^ is that it involves two spin- 
flip processes, created on adjacent copper-oxide planes. 
Since our present study relies on the Hubbard model on 
a single square-lattice plane, that mechanism would be 
beyond our theoretical approach. 

We recall that for each excitation momentum k, states 
in the (k, iS) domain of Fig. 2]whose energy is lower than 
that of the spin-wave spectrum intermediate sheet gen- 
erate no spectral weight and thus do not contribute to 
the spin dynamical structure factor. Furthermore, con- 
sistent with the first-moment sum rules of an isotropic 
antiferromagnet, no and nearly no weight is generated 
by states of any energy and momentum [0,0] = [0, 27r] = 
[2-7r,0] = [27r,27r] and near it, respectively. That to- 
gether with the small amount of spin spectral weight 
reported in Table [H] for energies between 450 meV and 
566 meV indicates that in that energy window there 
is nearly no spin spectral weight near the momentum 
[0, 0] = [0, 2tt] = [2tt, 0] = [2tt, 2tt] (see Fig. HJ). 

In addition to the Mott-Hubbard gap magnitude de- 
pendence on U /t, DMRG calculations were performed 
to derive the U/t dependence of the ground-state elec- 
tron single occupancy expectation value (1 — d). That 
quantity plays an important role in our study, in that 
it controls several spin-weight sum rules. Our prediction 
for the amount of total spin spectral weight in the energy 
range ui £ (0 meV, 450 meV) quantitatively agrees with 
that observed in the recent high-energy inelastic neutron 
scattering studies of Refi^, which were limited to that 
energy window. 

Moreover, as reported in Table HTl we predict that there 
is a small amount of extra weight ~ 0.1 \j? B above 450 
meV, which extends to about 566 meV. Since at and near 
the momentum [0, 0] = [0, 2tt] = [2tt, 0] = [2tt, 2tt] there 
is nearly no spin spectral weight, analysis of Fig. 5] re- 
veals that for energies between 450 meV and 566 meV 
the small amount of extra spin spectral weight is located 
at and around the momentum [tt, tt]. Thus we suggest 
that future LCO neutron scattering experiments scan the 
energies between 450 meV and 566 meV and momenta 
around [tt, tt]. 
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Appendix A: Useful operators algebra 

Here we justify why the six generators of the global 77- 
spin and spin SU(2) symmetries commute with the elec- 
tron - rotated-electron unitary operator. Moreover, we 
address the problem of the c fermion operator, spinon 
operator, and ?7-spinon operator algebras. 

To achieve our first goal, it is useful to express the 
kinetic-energy operator T given in Eq. ([1]) as T = To + 
f +1 +f_i. Here, 

T-y = - X] T T „j t j>, 7 = 0,±1, 
U,f) 

= ^"^fofj.-CT cj^.g- Cry, it ^ry ,— a 
a 

+ (1 ~ h^-a)ct (1 - fir'., _ CT ) +C.C.] , 

T+l;j,j> = ct. iff Cp.,^ (1 - hp.,- a ) 

a 

~l~ ^"r,/ ,— cr ., a c fj,a (1 ^"fj,— cr)] ; 



) C fj,a c fji,(r ^iy,— tr 



cr 

~f~ ( 1 ■/ ,— ct) ^V./ cr s<7 , — cr] 



(Al) 



While the operator To does not change electron double 
occupancy, the operators T+i and T_i change it by +1 
and —1, respectively. 

For U/t > the operator S in the expression V = e~ s 
given in Eq. (|43|1 can be expanded in a series of t/U, 



S 



t 

u 



0{t 2 /U 2 ) 



(A2) 



Although as discussed in Ref there are infinite choices 
for the operators V = e~ s and S, they share two im- 
portant propertie a 14 i 31 i 32 : (i) To leading order in t/U all 
read [f+i -T_i], as given in Eq. (|A2|) : (ii) Their op- 
erational expressions involve only the kinetic operators 
Tb, T+i, and T_i of Eq. (|A1 [) . Such properties apply to 
the specific electron - rotated-electron unitary operator 
V uniquely defined in this paper. 

The rotated kinetic operators To, J+i, and T_i such 



Tv = 



^0; 



7+1;; 



T-i, 



X! 7 " : • 7 = 0,±1, 

a 

(1 - _ CT )ct a gff , )tr (1 - + c.c/ 

^ ^ r,-, — cr ^Fj,cr -°" ~~ ^fy ,— cr) 
cr 

Tip., — n - C: 



r,,cr ^?j.cr (1 ™"Tj, — cr)] j 



X[(l ~ n^'-cr) ct Cp., <(7 hp., - a 



+ (1 - TV, -„) ct Cp^o-hp^-a] ■ 



(A3) 



To confirm that the three generators of the spin SU(2) 
symmetry, three generators of the 77-spin SU(2) symme- 
try, and also the momentum operator P commute with 
the electron - rotated-electron unitary operator V = V, 
one uses the exact result that the unitary operator V can 
be solely expressed in terms of the three kinetic opera- 
tors given in Eq. ^^ML. In RefJi the following 21 
commutators were found to vanish, 



[P,T 7 ] = [§Z,TJ = [Si,f 7 ] = [S a ,Tj =0, 

a = 77, s , 7 = 0, ±1 . (A4) 



Although the algebra involved in their derivation is cum- 
bersome, it is straightforward. The vanishing of the com- 
mutators given in Eq. (|A4[) then implies that the mo- 
mentum operator and the six generators of the 77-spin 
and spin algebras commute with the unitary operator V, 



[P,V] = [S z a ,V] = [SlV] = [S a ,V} = 0, 

a = 77, s, l = 0,±l. (A5) 



Hence the above operators have the same expression 
in terms of electron and rotated-electron creation and 
annihilation operators, so that the momentum operator 
reads, 



°-=t.4. g 



X 

<T=t4 k 



Furthermore, the above-mentioned six generators are 
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given by, 



I! / y a Tj,T) °Tj, f] 1 

3=1 J=l 

N a N a 

i=i 3'=i 

$V = X^,»7 = X] 5 ^,')' 

i=i i=i 

3=1 3 = 1 

3=1 3 = 1 

S ' = XV. Xv,- (AT) 

3=1 3 = 1 

Those of the unrotated local operators appearing here 
associated with the 77-spin algebra read, 



l r „ , 



c _ 



f ^c ft>t c r ~ M , j = l,2,...,N a , (A8) 



whereas those associated with the spin algebra are given 
in Eq. (j4j). On the other hand, the rotated local opera- 
tors appearing in the alternative expressions of Eq. (|A7[) 
are provided in Eq. 

The six local operators given in Eqs. (|3| and 
together with the local operator, 



<T=t4 



(!-%,_„), j = 1,2, ...,N a , (A9) 



are the seven generators of the U 7^ local gauge 
SU(2) x 5C/(2) x U(l) symmetry of the Hubbard model 
on a bipartite lattice with vanishing transfer integral, 

t = oia. 

Since the electron - rotated-electron transformation 
generated by the operator V is unitary, the rotated- 
electron operators ct. a and c?. >(T of Eq. (|43|) have the 
same anticommutation relations as the corresponding 
electron operators ct. a and Cf. jCT , respectively. Similarly, 

the local c fermion operators of Eq. ([50)1 and three lo- 
cal spinon operators and three local 77-spinon operators 
of Eq. have the same algebra as the corresponding 
unrotated spin-less and 77-spin-less fermion operators of 
Eq. and and three local spin operators of Eq. (|D) 

and three local 77-spin operators of Eq. (|A8J| , respectively. 
The former operators play a major role in the finite-///?; 
physics of the model. The latter operators are a limit- 
ing case of the former operators reached for U/t ^> 1 . 



Hence, without loss of generality in the following we pro- 
vide the algebra of the local c fermion operators of Eq. 
([50)1 and three 77s quasi-spin operators of Eq. ([S3)) . The 
SU (2) algebra of the latter three operators fully deter- 
mines those of the three local spinon operators and three 
local 77-spinon operators of Eq. ([S3]). 

Straightforward manipulations based on Eqs. (|50[) - 
(JSSJ lead to the following algebra for the c fermion oper- 
ators, 

{ff jt c > /ry,c} — j 

{4 ,c > 4' ,°> = {f? > ^ ' *v = ' (A10) 

and the c fermion operators and the local 77s quasi-spin 
operators, 



ft 

Jr. j 

ft si 



frj,c , Qr., 

f- si 

J rj 7 c ; "V -/ 



I = ±,X 3 . 



= 0, 



77, s. 



(All) 



The SU(2) algebra obeyed by the local 77s quasi-spin 
operators q\t. where I = £3, ±, such that <ji = g5. ± % q*L , 
and corresponding 77-spinon (a = 77) and spinon (a = s) 
operators si. a is, 



and 



2 <K 3 • 



T^j-' 9ft , (A12) 



St 



t* ■ rv " r ., 



°3,3 v a,a' * s F,-,a J 



a, a' = 77, s , (A13) 

respectively. Moreover, one has obviously that 
[9^ >«?.,] = and [«f!f,oi4.,,a'] = wnere i = 0,± and 
a, a' = 77, s. While the c fermion and 77s quasi-spin op- 
erator algebras refer to the whole Hilbert space, those 
of the 77-spinon and spinon operators correspond to well- 
defined subspaces spanned by states whose value of the 
number 25c of rotated-electron singly occupied sites is 
fixed. This ensures that the value of the corresponding 
77-spinon number M v = [N a — 2S C ] and spinon number 
M s = 2S C is fixed as well. 

The relations given in Eqs. (|A10p - (|A13[) confirm that 
when acting onto the model's Hilbert space the c fermions 
associated with the global c hidden U(l) symmetry are 
?7-spinless and spinless fermionic objects. They are con- 
sistent as well with the spinons and 77-spinons being spin- 
1/2 and 77-spin-l/2 objects, respectively, whose local op- 
erators obey the usual corresponding SU(2) algebras. 



Appendix B: Spin-wave spectrum in the high 
symmetry directions 

In this appendix we study the spin-wave spectrum of 
Eq. (|75p in the BZ high symmetry directions. These di- 
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rections correspond to those measured by high-resolution 
inelastic neutron scattering in LCO, as plotted for in- 
stance in Fig. 3 (A) of RefA (Our theoretical spin-wave 
spectrum curves are plotted along with the more recent 
LCO high-energy neutron scattering points of Refi£ in 
Fig. H) 

We denote such symmetry directions by MO, TO, 
XM, TX, and XO. They connect the momentum-space 
points M = [tt,tt], O = [7r/2,vr/2], T = [0,0], and 
X = [tt, 0] of the general spin- wave spectrum provided 
in Eq. ([75|) . The use of that equation reveals that the 
spin-wave excitation spectrum is in such symmetry di- 
rections given by, 



wro(fc) 
k 



sin(fc l ) , 



2 

[tt, -7r] - q - q" 

[h,h], ki = k x = k v £ (0,7r/2), (Bl) 



for si fermion hole momenta, 

q = [tt/2 - h, -tt/2 - h] , h G (0, tt/2) , 

= [tt/2, -tt/2], (B2) 



fe = [vr, 7r] — g — g 

= [fcj, fcj] , fc 4 = fc x = fey G (7r/2, 7r) , (B3) 

for si fermion hole momenta, 

<? = [tt/2 - fcj, 3tt/2 - h] , fei e (tt/2, tt) , 
g' = [tt/2, -tt/2], (B4) 



u)rx(k) 



sin(fc x /2) , 



(B5) 



fc = [tt, — tt] — q — q' 

= [*x,0], fc. G (0,7T), 

for si fermion hole momenta, 

q = [tt/2 - k x /2, -tt/2 - k x /2] , k x G (0, tt) , 
g" = [tt/2 - k*/2, -tt/2 + k x /2] , k x G (0, 7r) ,(B6) 



uxhi{k) = 
fc = 



cos(fc y /2) , 



[TT, fcj,] , k y G (0, TT) , 



for si fermion hole momenta, 



and, 



<f = [-fcy/2, TT - ky/2] , fey G (0, tt) . 

g' = [ky/2, -ky/2], ky € (0,7r), 



^xo(fc) = y - W^j cos(fc x ) 



+ W sl COS(ky) 



[tt, -7r\-q-q 
\TT,TT\-q" -q'" 



fox] i &x G (tt/2,tt) 
[tt - ky,k y ] , k y G (0,7r/2) , 



for si fermion hole momenta, 
q = [0, -tt] , 



[Q,TT 

Zy, ' 



[kn, k 



k y e (0,tt/2), 



(B7) 



(B8) 



(B9) 



[tt -k x ,-TT + k x ], k x G (7r/2, 7r) , (B10) 



(Bll) 



respectively. 

The theoretical spin excitation spectra, Eqs. (|Bip ~ 
(|B9|) . are plotted in Fig. [5] (solid line) for U/t m 6.1 and 
t w 0.295 eV together with the experimental results of 
Refi^ (circles) for T = 10 K. Such t//t and t magnitudes 
correspond to fj, — 565.6 meV and — 49.6 meV in 
the above energy spectra. The spin-spectrum expressions 
provided in Eqs. (fBT|) - (fB9|) refer to the first BZ. In Fig. 
[S]they are plotted in the second BZ, alike in Fig. 3 (A) 
of RefA An excellent quantitative agreement is reached 
for the above magnitudes of the involved energy scales. 
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